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Abstract. Consider the TJ^'^^-critical Schrodinger equation with a focusing cubic non- 
linearity in R"^ 

idt4^ + Alp +\ip\^ 4^ = 0. (0.1) 
It admits an eight-dimensional manifold of periodic solutions called solitons 

w(7r) = w{a, r,Vk,Dk)^ ^r(r+v.^-tH^+^H)^f^^ _ _ D,a), (0.2) 
where (f>{x, a) is a positive ground state solution of the semilinear elliptic equation 

-A(fi + a'^(fi^(fi^. (0.3) 

We prove that in the neighborhood of the soliton manifold there exists a H^^^ real- 
analytic manifold TV of asymptotically stable solutions of (IClf) . meaning they are the 
sum of a moving soliton and a dispersive term. 

Furthermore, we show that a solution starting on J\f remains on TV for all positive time 
and for some finite negative time and that A/" is a centre-stable manifold for this equation. 
The proof of the nonlinear results is based on the method of modulation, introduced by 
Soffer and Weinstein and adapted by Schlag to the L'^-supercritical case. Novel elements 
include a different linearization and a Strichartz estimate for the time-dependent linear 
Schrodinger equation. 

The main result depends on a spectral assumption concerning the absence of embedded 
eigenvalues in the continuous spectrum of the Hamiltonian. 



1. Introduction 

1.1. Main result. Consider equation (jO.ip : 

idfil; + Alp + = 0. 

For a parameter path tt = {v^, D^, a, T) such that ||'!T||iinL°° < ^^"^ ^ positive solution 
a) of dES]) 

define the nonuniformly moving soliton Wj^^t) by 

e{x, t) = v{t) ■x + (a^(s) - \v{s)\^) ds + T{t) 
y{t) = 2 [ v{s)ds + D{t). 
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Theorem 1.1 (Main result). There exists a codimension-one real analytic manifold J\f C 
H^l"^ , in a neighborhood of the soliton manifold, such that for any initial data ^(0) G M, 
equation i0.1\) has a global solution M and ip have the following properties: 

(1) ■0 is asymptotically stable, in the sense that it decomposes into a dispersive term 
r and a moving soliton WT^{t), given by 11. which converges to a final state: 

Ip = r + WTr{t) 

||^||i<Ca(0)||r(0)||^i/2 (1.2) 
||r|| -1/2 2Ti/i/2,6 < C'||r(0)|Ui/2 < C min ||'0(O) -w\\^i/2. 

-"a: I l^t ^'^^ -"a; W IS a Sohton -"a: 

(2) -0, r, and Wn depend real- analytically on the initial data ip{0)- 

(3) The dispersive term r scatters like the solution of the free Schrddinger equation: 
there exists tq € H'^^'^ such that 

r{t) = e''\Q + o^,/,{l). (1.3) 

(4) -0 stays on N for infinite positive time and for finite negative time, meaning that 
N is invariant under the Hamiltonian flow. 

(5) M is the centre-stable manifold of the equation. 

This result depends on the absence of embedded eigenvalues within the continuous 
spectrum of the linearized Hamiltonian (Section II. 5p . For a definition of centre-stable 
manifolds, the reader is referred to Section II. 3| for a definition of the norms involved in 
the statement of Theorem II. H one is referred to the Appendix. 

A similar result holds going backward in time. The invariant manifold in that case is 
obtained by conjugating the one described in Theorem ll.il 

As a final matter concerning notation, we denote by C various constants that appear 
in the proof, not all equal. 

1.2. Background and history of the problem. From a physical point of view, the 
nonlinear Schrodinger equation in with cubic nonlinearity and the focusing sign (jO.ip 
describes, to a first approximation, the self-focusing of optical beams due to the nonlinear 
increase of the refraction index. As such, the equation appeared for the first time in the 
physical literature in 1965, in [Kelley] . Equation (jO.ip can also serve as a simplified 
model for the Schrodinger map equation and it arises as a limiting case of the Hartree 
equation, the Gross-Pitaevskii equation, or in other physical contexts. 

More generally, consider the semilinear focusing Schrodinger equation in M"'^^ 

idtip + Ai; + lipl^ip = 0, ij{0) given. (1.4) 

It admits soliton solutions, e**"^(/)(x, a), where 

- A(/) + a^^ = (/>P+^ (1.5) 

Important invariant quantities for this equation include the mass 

M[iIj] = [ \i;{x,t)\'^ dx (1.6) 

and the energy 

Em= [ l\v^{x,t)\^-^mx,t)r'dx. (i.?) 
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Equation (jl.4p is invariant under the rescaling 

t) o^/Pilj{ax, Q^t). (1.8) 

We interpret this as meaning that ()1.4p is i/^'^-critical, for Sc = n/2 — 2 /p. Of particular 
interest are the Li^ or mass-critical (p = 4/n) and the or energy-critical {p = 4/(n — 2)) 
cases. 

Except for this introductory discussion of other authors' results, we always assume that 
n = 3 and p = 2, in which case (jl.4p reduces to (jO.ip and (jl.Sp to (j0.3p . In particular, 
(10. ip is i?^/^-critical; in general, p = 4/(n — 1) is i^^/^-critical. 

Given a soliton solution (II. 5p of the Schrodinger equation (II. 4p . a natural question 
concerns its stability under small perturbations. This issue has been addressed in the 
L^-subcritical case by Cazenave and Lions |CALlj and Weinstein |WEl2j . |WEl3j . Their 
work addressed the question of orbital stability and introduced the method of modulation, 
which also figured in every subsequent result. 

A first asymptotic stability result was obtained by Soffer-Weinstein [SoWElj . |SoWE2j . 
Further results belong to Pillet-Wayne [PlWAj . Buslaev-Perelman [BuFeI] . |BuPe2| . 
|BuPE3j . Cuccagn a [CuclJ, [ Cuc2], R odnianski-Schl ag-Sof fer [RoScSolj . |RoScSo2j . 
Tsai-Ya u [TsYAlj . |TsYA2j . |TsYA3j . Gang-Sigal [ GaSi] . and Cuccagna-Mizumachi 
ICuMlj . 

Grillakis, Shatah, and Strauss [GrShStT] , |GrShSt2] developed a general theory of 
soliton stability for Hamiltonian evolution equations, which, when applied to the Schrodinger 
equation, shows the dichotomy between the L^-subcritical and critical or supercritical 
cases. 

In the L^-supercritical, i?^-subcritical case in M^, Schlag proved the existence of a 
codimension-one Lipschitz manifold oiW^'^r\H^ initial data that generate asymptotically 
stable solutions to (jO.ip . This was followed by more results in the same vein, such as 
Buslaev-Perelman |BuPe1] , Krieger-Schlag |KrSc1] , Cuccagna [Cuc2j , Beceanu |Bec1| , 
and Marzuola [Mar] . 

In the L^-critical or supercritical case, negative energy {x)~^L'^nH^ initial data leads to 
solutions of (II. 4p that blow up in finite time, due to the virial identity (see Glassey [Gla] ). 
For a relaxation of this condition and a survey of results see |SuSu] and [Caz] . Berestycki- 
Cazenave [BeCa] showed that blow-up can occur for arbitrarily small perturbations of 
ground states p.Sp . Recent results concerning this subject include Merle-Raphael ^MeRa| 
and Krieger-Schlag [KRSc2j . 

Weinstein [WeiI] showed that all solutions (j) of the mass-critical equation ()1.4p of 
mass strictly less than that of the soliton, M[ip] < M[(j)], have global in time existence. 
Merle [MERj classified threshold solutions (M['0] = M[(^]), finding necessary and sufficient 
criteria for blow-up for data (namely, all blow-up solutions arise from transformations 
of the soliton), and global existence and scattering, for {x)~^L?' n data. 

A comparable result was obtained in 2006 by Kenig-Merle [KeMe] for the energy- 
critical equation ()1.4p in the radial case. Namely, for radial data </> of energy strictly less 
than that of the soliton, E[il)\ < £'[(/)], the following dichotomy takes place: if ||V'0(O)||2 < 
||V(/>||2, then (j) exists globally and scatters, while if ||VV'(0)||2 > \\^<P\\2 and V(0) G L"^ 
then ■0 blows up in finite time. In this regime, the equality ||VV'(0)||2 = ||V0||2 cannot 
occur. The behavior of solutions at the energy threshold, Elip] = E[(f)], was then classified 
by Duyckaerts-Merle |DuMe| : the same two cases are present, together with three more. 
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Following this approach, Holmer-Roudenko [HoRoj . Duyckaerts-Holmer-Roudenko 
[DuHoRo], and Duyckaerts-Roudenko |DuRo] established corresponding results for the 
ij^/^-critical equation (jO.ip . Their main findings may be summarized as follows: 

Theorem 1.2. Assume that ip is a solution of / (O.l)) with 

M[ij]E[ijj] - 2P[^f < M[<t)]E[^], (1.9) 

where (j) > is a soliton given by 10. 3\) and P is the momentum 

P[i;]= lp{x,t)iVij{x,t)dx. (1.10) 

Then one of the following holds: 

(1) //M[-0]||VV'||i-2P[V']^ < Af[0]||V(/)||2, thenip exists globally and scatters or equals 
a special solution, (p-, up to symmetries: Galilean coordinate changes, scaling, 
complex phase change, or conjugation. 

(2) If M[^p]\\V^\\l - 2P[^]^ = M[(p]\\V(p\\l, then ^ equals e'V(-,l) W to symmetries. 

(3) If M[ijj]\\Vilj\\l - 2P[tp]'^ > M[j)]\\V(j)\\l and V G {x)-'^L^ is radial, then blows 
up in finite time or must equal, up to symmetries, a special solution 

The special solutions 4>- and were defined by the aforementioned authors to have the 
following properties: (j)-{t) scatters as t ^ — cxd and converges at an exponential rate to a 
soliton solution as t ^ +oo, while p+it) blows up in finite time for t < and converges 
at an exponential rate to a soliton solution as t — > +oo. 

We explore the connection between their result, Theorem II. 2 1, and Theorem 1 1.1 1 in the 
subsequent remark. 

Remark 1.3. The boundary of the region 

B = {'4j\ M[iIj\E[iIj\ - 2P[ijf < M [</.]£;[</.]}, (1.11) 

is not a smooth manifold. Theorem 11.11 shows that, in the neighborhood of the soliton 
manifold, B is contained between two transverse hypersurfaces. The boundary of the region 
B and each of the two hypersurfaces intersect along nine- dimensional manifolds. 

The intersections are described by the special solutions 4'+{t) and (j)-{t), for sufficiently 
large t, subject to the symmetry transformations. 

This remark fits with several other natural observations. Firstly, eliminating the soliton 
manifold {(j) and its symmetry transformations — an eight-dimensional set) from B, an 
infinite-dimensional set, divides the latter into two disconnected components; this certainly 
could not happen if B were smooth. 

Secondly, consider the functional that defines B, 

F[7p] = M[7p]E[7p] - 2P[ipf. (1.12) 

Due to the extremizing property of the soliton cp in regard to the Gagliardo-Nirenberg- 
Sobolev inequality, the first differential of F at <j) is identically zero: 

lima*±itZM=0. (1,13) 

e^O e 

Therefore, the tangent cone to at (/> is actually given by the sign of second differential 
of F, which is an indefinite quadratic form; hence the lack of smoothness. 
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The most directly relevant results to which Theorem 11.11 should be compared are those 
of Schlag (SCHj , Beceanu |Bec1] , and Cuccagna j Cuc2] . 

In |Sch] . Schlag extended the method of modulation to the L^-supercritical case and 
proved that in the neighborhood of each soliton there exists a codimension-one Lipschitz 
submanifold of //-"^(R^) n VF^'"'^(M^) such that initial data on the submanifold lead to global 
n W^'^ solutions to (jO.ip . which decompose into a moving soliton and a dispersive 
term. 

[BecI] showed that for initial data in S = n H^, on a codimension one local 

Lipschitz manifold, there exists a global solution to (10. ip in the same space S. Further- 
more, the manifold is identified as the centre-stable manifold for the equation in the space 
S (in particular, the solution stays on the manifold for some positive finite time). 

Cuccagna [Cuc2] constructed asymptotically stable solutions for the one-dimensional 
mass-supercritical Schrodinger equation ( (11.40 for n = 1, 5 < p < oo) 

iut + Uxx + \u\'^u = 0, 5 < p < oo, (1-14) 

starting from even initial data (p = 5 is the L^-critical exponent in one dimension, 
while every exponent is iJ^-subcritical). The set of solutions obtained was not endowed 
with a manifold structure. 

Unlike previous results, the current one (Theorem ll.ip holds in a critical space, if 1/2, for 
equation (jO.ip . Previously it was not known whether the asymptotically stable manifold 
exists in the critical norm or is a phenomenon related to using stronger norms than the 
critical one in the study of the equation. The current work puts this question to rest. 

In addition, working in the critical space permits a series of improvements in the nature 
of results obtained. Firstly, we identify the asymptotically stable manifold as a centre- 
stable manifold for the equation in the sense of [BAJoj . Secondly, using the fact that the 
critical norm does not grow with time, we can prove that the manifold is globally in time 
invariant: solutions starting on the manifold exist globally and remain on the manifold, 
for all positive time. 

Finally, showing that the centre-stable manifold is real analytic raises the issue of its 
analytic continuation beyond the immediate neighborhood of the soliton manifold. It 
becomes an interesting problem to identify the global object that corresponds to the local 
centre-stable manifold. 

1.3. The centre-stable manifold. In 1983, Keller [Keller] considered the semilinear 
wave equation in with damping 

utt + aut- Au + f{u) = (1.15) 

and a stationary solution given by 

- Ano + /(no) = 0. (1.16) 

Under some growth and smoothness conditions on /, namely 

Fl lis) -ms = f{s) G C\R), m > 0; 

F2 /(O) = /'(O) = 0; 

F3 |/'(.)|<C(l + |s|2/(^-2)); 

F4 F{s) = f{t) dt<0 for some s > 0, 
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Keller proved the existence of an infinite-dimensional invariant local Lipschitz manifold 
of solutions that approach liQ as t ^ oo and of a finite-dimensional invariant local 
Lipschitz manifold of solutions that approach uq as t goes to — cxd. 

In 1989, Bates- Jones |BAJoj proved for a large class of semilinear equations that the 
space of solutions decomposes into an unstable and a centre-stable manifold. Their result 
is the following: consider a Banach space X and the semilinear equation 

ut = Au + fiu), (1.17) 

under the assumptions 

HI A : X ^ X is a closed, densely defined linear operator that generates a Cq group. 

H2 The spectrum of A decomposes into a (A) = as{A) U crc{A) U (Tu{A) situated in the 
left half-plane, on the imaginary axis, and in the right half-plane respectively and 
(7s{A) and (7u{A) are bounded. 

H3 The nonlinearity / is locally Lipschitz, /(O) = 0, and Ve > there exists a neigh- 
borhood of zero on which / has Lipschitz constant e. 

Moreover, let X", X^, and X^ be the ^-invariant subspaces corresponding to cr^j, cic, and 
respectively as and let S'^{t) be the evolution generated by A on X^. Bates and Jones 
further assume that 

Cl-2 dimX'', dimX" < oo. 
C3 Vp > 3M > such that ||S'^(t)|| < Me/'l*!. 

Let T be the flow associated to the nonlinear equation. We call J\f <Z U t-invariant if, 
whenever T{s)v G [/ for s G [0,t], T[s)v G for s € [0,t]. 

Let be the set of u for which T{t)u € U for all t < and decays exponentially as 
t — > — oo. Also, consider the natural direct sum projection tt'^^ on X'^ © X*. 

Definition 1.1. A centre-stable manifold N C U is a Lipschitz manifold with the property 
that M is t-invariant relative to U, 7r'^'*(AA) contains a neighborhood of in X"^ (B X^ , and 

MnW = {0}. 

The finding of |BaJo| is then 

Theorem 1.4. Under assumptions H1-H3 and C1-C3, there exists an open neighborhood 
U of zero such that is a Lipschitz manifold which is tangent to X'^ at and there 
exists a centre-stable manifold W^^ C U which is tangent to X^^ . 

Gesztesy, Jones, Latushkin, Stanislavova |GJLS] proved that Theorem 11.41 applies to 
the semilinear Schrodinger equation. More precisely, the authors established a spectral 
mapping theorem for the semigroup generated by linearizing the equation about the stand- 
ing wave. In particular, this theorem shows that the spectral condition 112 implies the 
semigroup norm estimate C3. The local Lipschitz property of the nonlinearity is not 
addressed in [GJLS| : it holds on for sufficiently large s. 

While providing a very general answer to the problem, the global existence of solutions 
on the centre manifold is not the subject of [BaJo| and [GJLS| . The i-invariance property 
of |BaJo] means that a solution starting on the manifold remains there for as long as it 
stays small. However, one does not know the global in time behavior of the solutions. 
Once a solution on the centre-stable manifold leaves the specified neighborhood of zero, 
not even its existence is guaranteed any longer. 
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By contrast, the estimates used by Keller [Kellejr] required a damping term, but his 
result held globally in time. Schlag |SCHj dispensed with the need for a damping term 
and proved a global asymptotic stability result, but the manifold that he constructed was 
not time- invariant. 

The current work identifies a centre-stable manifold for (jO.ip in the critical space for 
the equation H^/"^ and shows that solutions starting on the manifold exist globally and 
remain on the manifold for all time. 

1.4. Setting and notations. Equation (jO.ip admits periodic solutions e^^"^ (p{x, a), where 

(p = (j){x,a) = a(p{ax,l) (1-18) 
is a solution of the semilinear elliptic equation (10. 3p 

-A(j) + a^(t) = (f. 

In particular, we restrict our attention to positive Lp' solutions, called ground states. They 
are unique up to translation, radially symmetric, smooth, and exponentially decreasing. 
Their existence was proved by Berestycki-Lions in [BeLi| . who further showed that solu- 
tions are infinitely differentiable and exponentially decaying. Uniqueness was established 
by Coffman |COFj for the cubic and Kwong [Kwoj and McLeod-Serrin [McSe] for more 
general nonlinearities. 

Equation (jO.ip is invariant under Galilean coordinate transformations, rescaling, and 
changes of complex phase, which we collectively call symmetry transformations: 

0(t)(/(x, t)) = e*(^+^-^-*l'^l')a/(a(x - 2tv - D),ah). (1.19) 

If V'(i) is a solution to the equation then so is Q{t)tjj{t), with initial data given by g(0)'(/'(0). 

Applying these transformations to e'*(/)(") 1)) the result is a wider eight-parameter family 
of solutions to (jO.ip 

0(t)(e*V(^, 1)) = e^(^+^-^-*l'^l'+"'*)(/)(x - 2tv - D,a), (1.20) 

which we call solitons or standing waves. 

In the sequel we denote by a capital letter the column vector consisting of a complex- 
valued function (written in lowercase) and its conjugate, e.g. 

On this space we consider the real-valued dot product 

{F,G) = {f,g) + {7,g)= [ {f{x)g{x)+l{x)g{x))dx. (1.22) 

We leave it to the reader to check that our computations preserve this symmetry, in 
the sense that all column vectors that appear in the sequel comprise a function and its 

ff 

conjugate, meaning they are of the form F = \ -j 

We look for solutions to (jO.ip that get asymptotically close to the manifold of solitons. 
More precisely, we seek solutions of the form 

*=M'.(..t)+fl(..*)=(-;|^;;;)+(:(^;;;), (1.23) 
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where W-jt = ( ^ | is a moving soliton and R = ( !1 | is a small correction term that 
\w^J ^ \rj 

disperses like the solution of the free Schrodinger equation as t ^ +00. 

We parametrize the moving soliton 10-,^ by setting, as in (jl.ip . 

The quantities that appear in this formula, 

a{t), T{t), vit) = iviit),V2{t),V3it)), and D{t) = L»2(t), ^3 W), (1.25) 

are called modulation parameters and tt = (a, T, v, D) is called the parameter or modu- 
lation path. There are eight parameters in total, since v and D have three components 
each. 

Due to the nonlinear interaction between w.,^ and r, the modulation parameters are not 
constant in general; they are time-dependent. However, in the course of the proof they 
are not allowed to vary too much. A minimal condition, which we impose henceforth, is 
that 

d,t,Vk,DkeLl (1.26) 

and are small in this norm. We do not assume any stronger rate of decay. This implies 
that the modulation parameters converge as t — > 00, at no particular rate, and that their 
range is contained within arbitrarily small intervals. 

1.5. Spectrum of the Hamiltonian and the spectral assumption. The proof of 
Theorem 11.11 is based on a fixed point argument. Linearizing the equation around a 
moving soliton produces a time-dependent Hamiltonian of the following form, for a soliton 
w^o determined by a modulation path vr" as in ()1.24p : 

^ _ /^A + 2|u;,o(t)|2 {w^o{t))^ \ ^ 
^-"^'^"l -inlAt)? -A-2Mt)\^)- (1-27) 

We can always reduce the study of idt + 7i-,ro{t) to that of idt + 7i, where 

/A-l+2c/>2(.,l) X 

-<A^(-,1) -A + l-202(.,i);- (1-28) 

The spectrum of this Hamiltonian has been extensively studied. In Section 12.21 we give 
a summary of the known facts. In addition, though, we must make the following spectral 
assumption: 

Assumption A. The Hamiltonian 7i has no embedded eigenvalues within its continuous 
spectrum. 

Such assumptions are routinely made in the proof of asymptotic stability results, as, 
for example, in [BuPeI] . [Cuclj . |RoScSo2] . or |SCHj . 

1.6. Linear results. In the process of establishing the main nonlinear result. Theorem 
11.11 a crucial ingredient turns out to be a Strichartz-type estimate for the solution of the 
time-dependent linear Schrodinger equation 

idtZ - iv{t)VZ + A{t)a2.Z + nZ = F, Z{0) given. (1.29) 
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Here 

Wi and W2 are real- valued. 

For the purpose of establishing the nonlinear result Theorem II. H we may as well con- 
sider only potentials V that are smooth and exponentially decaying. 

We are mainly interested in Strichartz estimates 

\\z\\L^,L%nLrLi < ciwzmh + (1.31) 

Note that the Hamiltonian (jl.28p is nonself adjoint and time-dependent, leading to 
specific problems that are absent in the selfadjoint case. 

In the scalar selfadjoint setting, Keel-Tao |KeTa] proved endpoint Strichartz estimates 
for the free Schrodinger and wave equations and introduced a general method for obtaining 
endpoint estimates (based on real interpolation) that has been useful in all similar contexts. 

The issue of selfadjointness matters because one needs to reprove the usual dispersive 
estimates concerning the Schrodinger equation. They do not follow in the same manner as 
in the selfadjoint case, where, for example, the unitarity of the time evolution immediately 
implies the boundedness. 

In |SCHj . Schlag proved — > L°° dispersive estimates for the Schrodinger equation 
with a nonselfadjoint Hamiltonian, as well as non-endpoint Strichartz estimates. Erdogan 
and Schlag [ErSc) proved bounds for the evolution as well. In [BEClj . endpoint 
Strichartz estimates in the nonselfadjoint case were obtained following the method of Keel 
and Tao. Finally, Cuccagna and Mizumatchi [CuMi] obtained the boundedness of the 
wave operators, from which all of the above can be inferred as a simple consequence. 

In the linear setting, the most difficult to handle are terms of the form 

(ait) - a(oo))o-3Z and {v{t) - v{oo))VZ, (1.32) 

where Z is the solution. Instead of using Strichartz estimates to handle them, we make 
them part of the time-dependent Hamiltonian and thus avoid the issue altogether. 

Our main result in this context, which improves upon and is distinct from previous 
estimates such as |RoSc| . is the following: 

Theorem 1.5 (see Theorem 13. 9p . Consider equation ( fl.29|) . forH = Hq + V as in ( fl.30|) 

and a potential V not necessarily real-valued: 

idtZ -iv{t)VZ + A{t)a3Z + nZ = F, Z{0) given, 

M - ^ \,(Wi W2\ 
^ V -A + fij ^ \-W2 -WiJ- 

Assume that V = V1V2 Oind Vi, V2 are sufficiently smooth and decaying: 

{x)''d'^Vj G L°° (1.33) 

for every m, n < N and some sufficiently large N . Also assume that ||^||oo o,nd ||f ||oo 
are sufficiently small (in a manner that depends onV ) and (t('Ho) contains no exceptional 
values ofTC. Then 

ll^-^^llL-Hy^nLjTyi/^'^ - c{\\R{0)\\hi/2 + 11-^11^1^^2+^2^^1/2,6/5). (1.34) 
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The decay assumptions used here are not sharp. Pc is the projection on the continuous 
spectrum of 7i. More detailed statements of these results and remarks concerning them 
follow within the paper. For statements that hold under sharp or almost sharp conditions, 
see |Bec2| and |Bec3| . 

2. The Nonlinear Result 

2.1. Deriving the linearized equation. Substituting = w-,^ + r in the original equa- 
tion dnn]), 

idt{w^ + r) + /:\{w.„ + r) + {w.„ + r){w^ + rf =0. (2.1) 
In keeping with (jl.24p . the soliton w-,^ is described by 

w^{t) = e*(rW+/o(°' W-I^WP) ds+v(tyx) ^(^^ _ ^ g ^(5) - D{t),a{t)) . 
Note that WT^{t) depends on the values vr takes on [0, i]. Indeed, if we define 

w{-k) = w{a, r, V, D) = e^(''+^-^V(a; -D,a), (2.2) 
then Wj^it) 7^ w{'K{t))\ in fact, 

w^{t) =w(^{t),T{t) + {a^{s) - \v{s)\'^)ds,v{t),D{t) + 2 v{s)ds^. (2.3) 

Expanding the equation accordingly, note that 

dtWTT = (r + — v'^)drWTr + adaWjr + vd^w-,^ — (2v + D)dDWiT (2.4) 

and 

Au;^ = Ae^^(^'*V(a; - y{t),a{t)) + 2Ve^^(^'*) V(/)(a; - y(t), a{t))+ 

+ e*^(^'*)A0(x-y(t),a(t)) (2.5) 

= (a^ — v'^)Wt^ + 2iv\7WTr — IiUttI^U^tt- 

Here we used the following notation: given a soliton w = w{Tr) as in (j2.2p . drw, daW, 
doi^w, and dv^.w are the partial derivatives of w with respect to these parameters and span 
the tangent space to the soliton manifold at the point w: 

drw = iw, daW = daW, dD,.w = dx^w, dv^.w = ixkW. (2-6) 

We use these partial derivatives to define the differential dT^w of the map vr 1— > w = w{7r). 
The differential dj^w maps a vector = {Sa, 6T, 6vk, SD^) in the tangent space at vr to a 
vector in the tangent space at w: 

{d.,^w)67r = {dT)drw + {Sa)daW + {5Dk)dDf.w + {5vk)dv^w. (2.7) 

This expansion results in the cancellation of the main term involving w-,^ . The equation 
concerning r becomes 

idtr + Ar + i(9^u;^)7r + (|r|V + r^wJ^ + 2\r\^w-„ + 2r|w^p + rwl) = 0. (2.8) 

Here dT^w^^ is the differential evaluated at the point Wt^. Since as noted previously w-,^ / 
w[-k), (d^f^Tr)?!" is not the same as the time derivative of Wj^. 
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We regard equation (j2.8p as a Schrodinger equation for r. By conjugating, we obtain 
an equivalent equation for r: 

idfT — Ar + idj^WT^TT — (|r|V + (rj^u^^r + 2|r|^uJ^ + 2r|^i;7rp + r{wTr)'^) = 0. (2.9) 

Since both equations involve both r and r, it is most convenient to solve them together 
as a system or, rather, to see the pair of equations as just one equation concerning the 

Adopting this point of view, the main terms assemble into a time-dependent, nonselfad- 
joint two- by- two matrix Hamiltonian 

^ _^ -12:„p) ^ . V,,, (2.10) 

whereas the other terms are better treated as the homogenous right-hand side of the 
equation: 

' i{dT^WT^)'k + \z\'^z + z^IUtt + 2|zpi(;^' 



column vector R 



dT^WT^TT — \z\'^z — {z)'^Wt^ — 2\z\'^WTrJ ' C^-^^) 

Given a C^-valued soliton W = (^^^ ' ^^^o introduce the partial derivatives djW = 

(§^) ' ^""^ ^ "l") r, Vk, Dk}: 

drW = (^^^^ = ia3W, daW = daW, do.W = d,,W, d,,W = ia^x^W. (2.12) 
As a reminder, is one of the Pauli matrices: 

-3=(; \). (2.13) 

We use these to define the differential dT^W^-K) of the map 

^ ' \w{tt)) (2.14) 
(d^Ty(7r))57r = {5T)dTW + {5a)dc,W + {5Dk)dD,W + {5vk)d,,W. 
Thus, the first half of F is simply given by the differential 

'i{d^w^)T:\ _ ^2.15) 



i dj^WT^'k 

We denote the remaining nonlinear term by 

N{R, W^) = ( - f - . (2.16) 

In vector form, the equation fulfilled by R becomes 

idtR - n^{t)R = F{t), F = -i{d^W^)n + N{R, W^). (2.17) 

To this equation concerning R we join the modulation equations that determine the 
path vr. For future convenience, define the cotangent vectors 

E^{W) = iasdrW, Er {W) = ia3daW, 

E,^iW) = iasdo.W, Ed,{W) = iasd^.W. ^^'^^^ 
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Also consider the (second) differential d.„'Ef{W), naturally defined as the differential of 
the map tt i-^ EfiW^n)), for / G {a,T,Vk,Dk}. 

At each time t and for all / G {a,T,Vk, D^} we impose the orthogonality condition 

{R{t),Ef{W^{t))) =0. (2.19) 

This condition arises as follows: given parameters vr = {a,T,Vk, D^), let 

w = w{Tr) = e'^^+''-''U{x - D,a), 

W = W{.) = Q , ^2.20) 

niw)^[^^_y _^:;,^p)+2..v-(a^-HV3. 

The Hamiltonian 7i{W), naturally associated to each soliton W (see Section l2.2p . will 
play an important role in the sequel. 

Note that n{W^{t)) = n^{t) + 2iv{t)y - {a{tf - |f (t)P)(T3. 

Condition (j2.19p is then equivalent to asking that R should not live in the zero eigenspace 
of the time-dependent Hamiltonian 'H(W^(t)). 

Taking the derivative in (j2.19p . it translates into the following modulation equations. 

Lemma 2.1 (The modulation equations). 

Q = Aa\\W^\\-^^{{R, {d^E^{W^))iT) - i{N{R, W^), E^{W^))) 
f = 4a\\W^\\^\{R, {d^EriW^))^) - i{N{R, W^),Er{W^))) 
Vk = 2\\W^\\f{{R, {d^E,^{W^))ir) - i{NiR, W^),E,^iW^))) ^^'^^^ 
Dk = 2\\W^\\^\{R, id^ED,{W^m) - i{N{R, W^),EDMn)))- 
Proof. Begin by observing that for every soliton W 



{daW, 


Ef{W)) = 






= a and zero otherwise 


{drW, 


Ef{W)) = 


4a 


\l if/ 


= r and zero otherwise 




Ef{W)) = 


2" "2 


if / = 


Dfc and zero otherwise 




Ef{W)) = 


2" "2 


if/ = 


Vk and zero otherwise. 



(2.22) 



Furthermore, 



Finally, 



a ) 



n*iW)Ea,iW) = 0, n*{W)Er{W) = -2iE, 
n*{W)E,^iW) = 0, n*iW)ED,m = -2iE,^. 



(2.23) 



dtEfiw^) = (d^Efiw^))^ + {n*{w^) - ni)Ef{w^). (2.24) 

Then, in the equality 

{R,dtEf{W^)) = -{dtR,Ef), (2.25) 

we replace dtR by its expression ()2.17p and arrive at (I2.2ip . □ 
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Let 



L^R=Aa \\W^\\^^{R,{d^Ef{W^))ir)dfW„ 

/Ga,r 

+2 Yl \\W^\\2^{R,{d^Ef{W^))n)dfW^ 

/6K-,Dfc} 



(2.26) 



and 



(2.27) 



N^iR,W^)=^a Y \\W^\\2^i{N{R,W^),Ef{W^))dfW^ 

fe{a,T} 

+2 Yl \\W^\\2^i{N{R,W^),Ef{W^))dfW^. 

The modulation equations can then be rewritten as 

{d^W^)n = L^R- iN^{R,W^). (2.28) 

L-^R represents the part that is hnear in R and Nt^{R,W) represents the nonhnear com- 
ponent {N{R,W),Ef{W^)). 

As an aside, using the notation of Section [2.21 let Po(^7r) be the zero spectrum pro- 
jection of the Hamiltonian corresponding to Wj^. We recognize 

W^) = Po{W^{t))N{R, W^) (2.29) 

and a similar, but more complicated, expression for L-,^R. 

At this point we linearize the equation. We use an auxiliary function for all quadratic 
and cubic terms and do the same in regard to the soliton. We only keep linear, first-order 
terms in the unknowns R and vr for which we solve the equation, so that the equation 
becomes linear in R and vr and quadratic and cubic in the terms involving R^ and vr^. 

We introduce an auxiliary function R^ and an auxiliary modulation path vr*^. Starting 
with vr*^, we construct the moving soliton w^o, its partial derivatives dfw^o, cotangent 
vectors EjiW^o)^ the differentials dT^W^o and d^Hj(VKn-o), etc., all following previous def- 
initions: 

w^o = e^(r"W+/o W-l-''PW)'^-+-''W-)</,(x -2j^v^{s)ds - (t) , (t)) , 



(2.30) 



as well as 



dfW^o 



dfW^o 



yw^oj (2.31) 



'^aiW^o) = icj3(9rVF^o,Hr(W^o) = iasd^W^o, 
Lemma 2.2. ^ is a solution of liO.l]) if and only if 
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W^o is given by 112.30]) . and {R,it) is a fixed point of the map {R^,tt^) ^ {R,7r), where R 
and TT solve the following system of linear equations: 

idtR + n^o{t)R = F, F = -iL^oR + N{R'^,WA - N^o{R°,WA 

f = 4aO||M^,o||2 '((i?, K%(T^,o))7^0) _ .^^(^0^ WA,Ef{WA)), f E {a, T} (2.33) 

/ = 2\\WA\^'{{R, {d^Ef{WA)n') - i{N{R^, WA,Ef{WA)), f G {vk,Dk}. 



Here 

'A + 2\wAt)? {wAt)? 

-{wAt)? -A-2|«;,o(t)|2'' 

- (r^w^o - 2\r^\^wA 

inn ,W^o)- I |^0|2^o + (^0)2 +2|r0|%^o 

I I TT , (2 34) 

L^ofl = 4Q° \\Wno\\2^{R,{d7rEf{WA)^^)dfW^o 
/6Q,r 

+ 2 Yl Il^-«ll2 id^EfiWA)^'')dfW^o, 

fe{vk,Dk} 

/6{a,r} 

+ 2 Yl \\WA\2MNiR^,WA,^fiWA)dfW^o. 

Solving the system (|2.33p in a suitable space, with a view toward applying a fixed point 
theorem to conclude that {R^,tt^) = (i?, vr), is our main objective. 

Proof. We substitute R^ for R and vr'^ for vr into all higher order terms of (|2.17p and 
()2.2ip . Following this substitution, the nonlinear Schrodinger equation (|2.17p transforms 
into the linear equation, in both R and tt, 

idtR-n^o{t)R = F, (2.35) 

where 

F = -i{d^WA^ + NiR'^, WA, 

(2.36) 

{d^WA^^ = TdrW^o + adaW^o + DkdD,,W^o + Vkdy^W^o. 
The orthogonality condition we impose in the linear setting is 

{R{t),Ef{WAm =0, fe{a,r,Vk,Dk}. (2.37) 

By taking the derivative, in a manner entirely similar to that of Lemma |2. II we obtain 
the linearized modulation equations: (j2.2ip become 

a = 4a''\\WA\2\{R, {d^Ea{WA)A " ^(^^(«°, vfO), 5° )) 

t = 4a^W%\{R, {d^Er{WA)A - i{N{R',7r^):E^r)) 

Vk = 2\\W%\{R,{d^EAWA)A-^{N{R',A,^V) 

Dk = 2\\wX\{R,{d^EDdWA)^')-i{NiR',A,^DA- 
The identity (j2.28p translates into 

{d^WA-^ = L^oR - iiV^o(i?°,7r°). (2.39) 



(2.38) 
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Finally, we collect together (j2.35p and (j2.38p and replace (dTrW^vr")'^ on the right-hand 
side of (j2.35p by its expression (j2.28p in order to arrive at (|2.33p . □ 



2.2. Spectral considerations. Given parameters vr = (a,T,vi:, Dk), define the soliton 

w = w{7r) = e'^'^+''-''Uix - D, a), 

W = W{7r) = 



(2.40) 



and consider the associated Hamiltonian 

n{W) = n{a, r, v, D) = _^ ^ ^l^^) + 2^W - (a^ - \v\'')a,. (2.41) 

By rescaling and conjugating by as well as by a translation, one sees that all 

these operators are in fact conjugate, up to a constant factor of a^: 

Dili/„e-^^TDe-^(^^+r)'^37^(a,r,?;,Z))e^(^^+r)'^3e^^Dil« = a27^(l, 0, 0, 0). (2.42) 

Therefore all have the same spectrum up to dilation and have similar spectral properties; 
thus, it suffices to study TL = TC{1, 0, 0, 0): 

/A-l + 2<^2 ^2 X 

_A + l-202j- (2-43) 

We restate the known facts about the spectrum of 7i. As proved by Buslaev, Perel- 
man |BuPe1) and also Rodnianski, Schlag, Soffer in [RoScSo2] . under fairly general 
assumptions, ct('H) C MUiM and is symmetric with respect to the coordinate axes and all 
eigenvalues are simple with the possible exception of 0. Furthermore, by Weyl's criterion 

(Tessin) = (-00, -1] U [1, +00)^ 

Grillakis, Shatah, Strauss [GrShStT] and Schlag [SchJ showed that there is only one 
pair of conjugate imaginary eigenvalues iticr and that the corresponding eigenvectors de- 
cay exponentially. For the decay see Hundertmark, Lee |HuLe] . The pair of conjugate 
imaginary eigenvalues ±ia reflects the L^-supercritical nature of the problem. 

The generalized eigenspace at arises due to the symmetries of the equation, which is 
invariant under Galilean coordinate changes, phase changes, and scaling. It is relatively 
easy to see that each of these symmetries gives rise to a generalized eigenvalue of the 
Hamiltonian Tl at 0, but proving the converse is much harder and was done by Weinstein 
in |Wei2], [Wei3] . 

Schlag |Sch| showed, using ideas of Perelman [Per) , that if the operators 

L± = -A + 1-2(I)\;1)T^\;1) (2.44) 
A i \ 

that arise by conjugating 7i with I . I have no eigenvalue in (0, 1] and no resonance 



at 1, then the real discrete spectrum of 7i is {0} and the edges ±1 are neither eigenvalues 
nor resonances. A work of Demanet, Schlag [DeSc] verified numerically that the scalar 
operators meet these conditions. Therefore, there are no eigenvalues in [—1, 1] and ±1 are 
neither eigenvalues nor resonances for TC. 

Furthermore, the method of Agmon |Agm| . adapted to the matrix case, enabled Erdogan- 
Schlag [ErSc] and independently [CuPeVo] to prove that any resonances embedded in 
the interior of the essential spectrum (that is, in (— oo, — 1)U(1, oo)) have to be eigenvalues, 
under very general assumptions. 
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Under the spectral Assumption [A] we now have a complete description of the spectrum 
of 7i. It consists of a pair of conjugate purely imaginary eigenvalues zizia, a generalized 
eigenspace at 0, and the essential spectrum (—00, —1] U [1, 00). 

Following |SCHj . let be the eigenfunctions of Tl corresponding to iticj; then there 
exists £ LP' such that 

Due to the symmetry 

asWda = n* (2.46) 

of the operator, the respective eigenfunctions of Ti* are a^F^ . Then, the imaginary 
spectrum projection is given by 

Pi^ = P+ + P_, P^ = (.,m3FT)F± (2.47) 

up to a constant; the constant becomes 1 after the normalization 

/ Re/+(x)Im/+(x)dx = -1/2. (2.48) 

It helps in the proof to exhibit the discrete eigenspaces ofn{W). Denote by F^{W) 
the eigenfunctions of W(VF) corresponding to the ±ia{W) eigenvalues. Even though there 
is no explicit form of the imaginary eigenvectors, Schlag |SCHj proved in a more general 
setting that F^{W), normalized, and a are locally Lipschitz continuous as a function 
of W and that F^ (W) are exponentially decaying. 

More is true in the case under consideration. Since the operators Ti.{W) are conjugate 
up to a constant, the dependence of F^{W) and a on the parameters can be made explicit: 

F^(W) = e*(^^+^)'^3e^^DiUF±, 

, (2.49) 
a{W) = a^a. 

Also observe that dfW^ where / G {a, r,r;fc,Dfc}, are the generalized eigenfunctions of 
n{W) at zero and S/(I^), defined as in (f2J[8]l . fulfill the same role for n{W)* . 

We then express the Riesz projections corresponding to the three components of the 
spectrum of TLiW) as 

Pim{W) = P+{W) + P^{W), P±{W) = a-^-MzF^{W))F^{W), (2.50) 

Po{W) = 4a||ir||2 '((•, '^a)daW + (•, ■Br)drW)+ (2.51) 



and 



Pc{W) = l-Pirr,{W)-Po{W). (2.52) 



2.3. The fixed point argument: stability. Consider a small neighborhood of a fixed 
soliton wq = w{'Kq), determined by parameters ttq = (aoi Tq, ffcOi ^ifco): 



.«^oy (2.53) 
^„o = e^(ro+-o-)</.(x-A,ao). 
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As in Section \2.2\ Wq has an associated Hamiltonian 7i(Wo) of the form (j2.4ip 

n{Wo) = +4^1' _^ + 2^voV - {al - \vo\')as. (2.54) 

In turn, to 7i(Wo) we associate the zero spectrum projection Pq{Wo), the imaginary spec- 
trum projection PimiWo) = P+{Wo) + P-{Wo), and the continuous spectrum projection 

PciWo) =1- PoiWo) - P^„^iWo). 

Up to quadratic corrections, the centre-stable submanifold is given by the affine sub- 
space 

Wo + {PciWo) + P.{Wo))H^/^ = {Wo + Ro\Ro£ H^'\ (Po(^o) + P+{Wo))Ro = 0}. 

(2.55) 

This submanifold will have codimension nine, so we need a supplementary argument, 
presented at the end, to recover eight codimensions. 
Take initial data of the form 

i?(0) = Ro + /iF+(Tyo), ^(0) = vro = ("o, Tq, {vk)o. {Dk)o), (2.56) 

where Ro £ {Pc + P-)H^^^ . Thus, Ro belongs to the appropriate linear subspace, while 
hF^ is the quadratic correction, made in the direction of F^; in particular, 

Po(0)i?(0) = 0, P+R{Q) = h. (2.57) 

(|236D also implies that ||i?(0)||^i/2 < C{\\Ro\\hi/2 + \h\). 
We define a map T as follows: 

Definition 2.1. T is the map that, given a pair {R^^-k^), associates to it the unique 
hounded solution (R, vr) of the linearized equation system i2.33\} from Lemma 12.21 with 
initial data as in i f2.56|) . where Rq is given and h is allowed to take any value: 

T((i?o,7rO)) = (i?,vr). (2.58) 

In the sequel we show that, for given Rq, R^, and n^, the bounded solution {R, vr) exists 
and is unique and that the parameter h = h{Ro, R^ ,tt^) is determined by the condition 
that the solution should have finite X seminorm. There exist a unique value of h for 
which the solution is bounded and a unique solution (R, vr) corresponding to that value of 
h. Thus the map T is well-defined. 

Consider the space 

X = {{R,7t) \ R e Lrn'J^ n L^^W^/^''', irGl'}, (2.59) 

with the seminorm 

||(i?,7r)||x = \\P\\looh1/2^^^i/2,6 -Fao^||7i-||^i, (2.60) 

where qq is the scaling parameter, see (j2.56p . X is the natural space for the study of 
equation (jO.ip and of its linearized version (j2.33p . Indeed, since the Schrodinger equations 
(10. ip and (I2.17P are if^/^-critical, we need to study them in the critical Strichartz space 
Lf^Hy^ n L^Wx^"^'^ ■ In regard to the modulation path, tt G is a minimal assumption 
to ensure we are dealing with small perturbations — and there is no room for a stronger 
condition, due to working in a critical space. 

We prove that, for \\{R^ ,tt°)\\x < 5, it follows that ||(i?,^)||x = \\T{R° ,tt^)\\x < 5 as 
well, when 6 is small. For clarity, we state the stability result formally: 
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Proposition 2.3. Let {R,tt) be a solution of < [2.33|) with initial data given by < [2.56p . 
Assume that the auxiliary functions {R^,tt^) of 12 .33)) also satisfy {2.56)) and, for some 
sufficiently small 5 < 5o, \\{R^ ,tt^)\\x < 5- Then (i?, vr) G X for a unique value of h, 
which we denote h{Ro, R^,tt^), and this bounded solution satisfies 

ll^llLr^^^'nL?H/iA6 < C(II^0|Ihi/2 + <5^), 
11^11^1 + |/i(i?o,i?°,vrO)| <CaoiS\\Ro\\Hi/2+S^). 
Supposing that ||i?o|liji/2 < c6 and 5 is small, it follows that the set 

{{R,7T)\R{0) = Ro + hF+{Wo), ^(0) = ttq, \\R,7t\\x<S} (2.62) 
is stable under the action of the map T and that 

\h{Ro, R°, Tr^)\ <C5\ (2.63) 

To a first order, R itself is then given by the time evolution of Rq under 7i^o{t) and all 
other terms in its composition are of higher order (of size 5^). 

We claim the previous result for each Wq and sufficiently small 6, but are also interested 
in the dependence of 6 on Wq. Since the Strichartz norms are scaling-invariant, they are 
left unchanged by symmetry transformations; it and h scale like ao- After accounting for 
this fact, we are left with an invariant statement. 

Proof. The modulation path vr'^ : [0, oo) M.^ determines the moving soliton w^o, follow- 
ing (j2.30p . and the time-dependent Hamiltonian 'H^o(t), as per Lemma l2.2[ 

Having in view Theorem 13.91 and the preceding discussion in Section 13.31 recall the 
notations 

w{7r\t)) = e*(r°(t)+^'°W-)^(^ _ D'{t),a''{t)), 

and consider the family of isometries 

^ Jo{2v'{s)V + .((a0(.))2 - \v'{sr)as) ds _ ^2.65) 

Observe that 

W{7:\t)) = U{t)WAt), 

n{w{7:\t))) = u{t)-^n{wAm{t), 

and let 

Z{t) = U{t)R{t). (2.67) 

For brevity we denote 



(2.66) 



Hit) 


= n{W{7r\t))), 


Po{t) 


= Po(TV(vrO(t))), 


Pc{t) 


= P,{W{n\t))), 


Pim (0 


= PUW{n\t))) 


P+it) 


= P+iWiTT^t))), 


P-{t) 


= P_(mvr°(t))), 


F+{t) 




F-{t) 


= F~{W{^\t))l 


Ef{t) 




a{t) 


= a{W{7r'm. 



(2.68) 



Rewriting the equation ()2.33p from Lemma 12.21 

idtR - n^o {t)R = F{t) (2.69) 
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with Z as the unknown, we obtain 

idtZ -n{t)Z = U{t)F{t). (2.70) 

Instead of fixing a time-independent Hamiltonian for this equation, we consider the 
time-dependent Hamiltonian Ti.{t) and divide the equation for Z into three parts, according 
to the three components of the spectrum — continuous, null, and imaginary. Let 

I = Pcit) + Poit) + Pi^it), Pi^{t) =P+(i)+P_(t). (2.71) 

Then, we separately prove, for each of the three components, estimates that enable us to 
carry out the contraction scheme. 

The Pq component is the most straightforward. By (j2.56p . the orthogonality condition 

{Rit),EfiWAm=0 (2.72) 

holds at time t = and the modulation equations (I2.33P then imply that the orthogonality 
condition still holds at any other time t. Applying the isometry U{t), the orthogonality 
condition turns into 

(Z(t),H^(t)) = 0, (2.73) 

which directly implies that 

Po{t)Z{t) = 0. (2.74) 
The continuous spectrum projection of Z fulfills the equation, derived from ()2.70p . 

idt{Pcit)Z) + n{t)Pc{t)Z = P,{t)U{t)F{t) + i{dtPc{t))Z. (2.75) 

The right-hand side term 

F = -iL^oR + N{R^, TT^) - N^o {R°, tt^) (2.76) 

is bounded in the dual Strichartz norm by means of the fractional Leibniz rule: 



l^7r0^llr2Ti/l/2.6/5 ^ ^^11^11^2^1/2,6 = C6\\Z 



r2wl/2,6, 



||iV(i?o,7rO) - iV,o(flO,7rO)||^2^v2,6/5 < C6\ ^^'^'^^ 

When using the fractional Leibniz rule, there is an endpoint Sobolev embedding issue, 
namely that W^^'^'^ does not embed into L°°. However, when estimating the quadratic 
and cubic terms, such as Ir^^r", present in N{R^, n^), we can avoid the problem by putting 
in non-endpoint Strichartz spaces: 

||(rO)2^(vrO)||^,^V2,6/5 < C||(r0)2|| ,/2||u;(vr0)||^^^v2 

t ^ t ^ t ^ /r\ '-7C)\ 

<C||i?l^,^./2,3ir(vr'')||,.^y2. 

The same applies to the cubic term, except that there we replace wlir^) by r^. 

The endpoint Strichartz space intervenes only when evaluating linear terms such as 
L^oR. However, there we avoid the attendant endpoint Sobolev embedding issue by means 
of a less sharp fractional Leibniz rule, which we can use because solitons are of Schwartz 
class. 
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Furthermore, 

{dtP+{t))Z = {dta'^){Z,iazF~{t))F+{t) + a-^{Z,icj^dtF-{t))F+{t)+ 
+ a-'{Z,ia3F-it))dtF+{t) 

= -3a-'^d{Z,ia3F~{t))F+{t) + a-^{Z,ia3d^F-{t)n^{t))F+{t)+ 
+ a-^{Z, ia3F~{t))d^F+{t)'fr'^{t) 
and likewise for P and Pq. Since Pc = I — Pq — Pim, by direct examination it follows that 

ll(5.P±(t))^||,.^./.a/s < C^'^ll^ll^.^iAa. (2.80) 

Provided ||vr''||i < 5 is sufficiently small, Theorem 1 3 . 9 1 leads to Strichartz estimates for 
Pc{t)Z, by means of the following construction. 'H(t) is given by 



-w'in^it)) -A-2|u>(7r"(t))|^ 

(2.81) 

Denote, for 7r''(0) = ttq following ()2.56p . 

The difference Ti — TC is small in the appropriate W^^"^'^^^'^ Pi norm (in much 

stronger norms too), so the corresponding term can be bounded by means of endpoint 
Strichartz estimates: 

\\{H - H)Z\\ ^2^^1/2,6/5 < \\T~C{t) - W(i)||^oo^l/2,6/5-e(-l^l/2, 6/5+^11-^11^2^^^2,6 

(2.83) 



< C||7r°(t) -7ro||L-||Z||^2^i/2,6 

< C|l^°llill^lli2^^i/2.6 

< C6\\Z\\ ^2^1/2,6. 

Then, by ()2.77p . (j2.80p . and (j2.83p . for sufficiently small 6 and under the assumption 
that there are no embedded eigenvalues we can apply Theorem 13.91 with the Hamiltonian 
Ti. and obtain 

\\Pcit)Z\\^^^ij2^^2^ij2,e<C{\\Zmm^^ (2.84) 
For the imaginary spectrum projection, write, following Section [2 ■2[ 

Pimit)Z{t) = b+{t)F+{t) + b'{t)F-{t). (2.85) 
Then, by taking the time derivative of 

6±(t) = (aO(t))-3(Z(t),icT3FT(t)), (2.86) 

we obtain 

dtb^ = {a'')~\dtZ,ia3F^) - 3a''{a'^)-\Z,ia3F^{t)) + 
+ {a'')-^{Z,ia3{d^F^)^'). 
Substituting dtZ by its expression given by equation (j2.33p . we arrive at 
= ±cT(i)6^ - {UF,a3F^)- 

- 3d°(a°)-^(Z,m3FT) + (a°)-3(Z,m3(d^FT)^°). 



(2.87) 



(2.88) 
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Thus and b satisfy the equation 
where 

N± = -{UF,a3F'^) -3cP{a^r'^{Z,ia3F^) + {a^y^ {Z,ia3{d^F^)n^) . (2.90) 

Here zizia{t) are the imaginary eigenvalues of TCi^t), as in our discussion of its spectrum in 
Section Eia 

Concerning the right-hand side, (j2.88p and (j2.77p imply 

\\N±{t)h2 < C{6\\Z{t)\\^,^,,,^e+6'). (2.91) 

To control the solution we use the following simple fact, see [Sch]. It characterizes the 
bounded solution of the ordinary differential equation system (I2.89j) . 

Lemma 2.4. Consider the equation 

'"o -°(())^- = /("' P'^^) 

where x = ^^^^ unknown, f = (jjf^ ^ ^ ^T' > '^'^'^ ^{t) > ctq > is bounded 

from below. Then x is bounded on [0, oo) if and only if 

= a;i(0)+ / e~^o^(^)<i^f^{t)dt. (2.93) 
Jo 

In this case, for all t >0 

xi{t) = - / efs-(-)^-f,{s)ds, 



. (2.94) 

a;2(t) = e-io-W'='-x2(0)+ / e-fs''(^)'i^f2(s)ds. 

Jo 

Proof. Any solution will be a linear combination of the exponentially increasing and the 
exponentially decaying ones and we want to make sure that the exponentially increasing 
one is absent. It is always true that 



XI (t) = e/o (xi (0) + e- f^ (s) ds) , 

(2.95) 

X2(t) = e-/o-W'^-X2(0)+ / e-/>(^)'^V2(s)(is. 



Thus, if xi is to remain bounded, the expression between parantheses must converge to 0, 
hence ST^ . Conversely, if dlTMl) holds, then 



oo 



xi{t) = -j^ e^s^''^)'i^f^(s)ds (2.96) 
tends to 0. □ 
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Note that a depends Lipschitz continuously on the scahng parameter a, in an expUcit 
manner, fohowing (j2.49p . Then a{t) belongs to a compact subset [ai, 02] of (0, 00), because 
a^{t) belongs to a compact subset of (0, 00). Consequently, equation (|2.89p has a bounded 
solution if and only if 

= 6+(0)+ / e-^o^M'^^iV+(t)dt. (2.97) 







One sees that 6+(0) = h, where 6+(0) is given by ()2.97p and h by (j2.56p . R is globally 
bounded in time if and only if Z is bounded. Z is bounded if and only if each of its 
components is bounded, PimZ in particular. Thus R is bounded only if 

/■CO ^ 

h = - e-^o^(^)'i^N+{t))dt. (2.98) 



Jo 

We are interested in a more direct formula for h, one that involves R instead of Z. Note 
that ct^t) is also the imaginary eigenvalue of TC(W^o{t)). Expanding A^-|_ and reverting 
the isometry U within (j2.98p leads to the explicit formula 

h = - / e-io"(^-°("»'^"((F,cj3F-(VF'^o(t)))- 
Jo 

-3a^{t)ia'it))~^{R,ia3F-iW^om+ ^^'^^^ 
+ {a\t))-^{R,ia-i{d^F-{WAt)))'fr'^{t))) dt. 

It remains to show that, for this unique value of /i, R is indeed bounded. 
We obtain 



\h\<C e~^o'^W'^^|iV_^(t)|dt 

<C e~'-^\N+{t)\dt (2.100) 
<C||iV+||ii+^c. 

Following (I2.94p . both 5+ and 6_ are given by convolutions with exponentially decaying 
kernels in t, whose rate of decay is bounded from below: 



|6+(t)| < / e^'-'^'^^\N+{s)\ds, 

).t (2.101) 

|6_(t)|< / e-(*-^)'^MiV_(s)|(is + e-*'^MI^o|liji/2, 
^—00 

with the convention that N{s) = for s < 0; the second term for h-{t) stems from 
e-*''6_(0). One has 

\\Pim{t)Z\\ ^2^i/2fi < \\h+\\Ll + 

< C{\\N+\\l2 + \\N_\\^2 + ||i?o||^i/2) (2.102) 
<C{\\R4hU^+5\\R\\^,^,j,^, + 5^). 

Putting all three estimates (I2.74j) . (I2.84j) . and (|2.102p together and taking into account 
the fact that 
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we obtain 

11^^11^2^1/2,6 < C{\\Ro\\hV2 + 5||ii||^2^1/2,6 + 5^). (2.104) 

Concerning the modulation path vr, from the modulation equations (j2.33p we get that 

\\ii\\i<C{6\\Z\\^,^,j,,,+6^). (2.105) 

Overall, 

\\{R,^)\\x = \\{Z,7,)\\x<C{\\R4^,j,+6\\{R,7i)\\x + 5^) (2.106) 
and this proves stability for the suitable choice of h and small initial data Rq. □ 

2.4. The fixed point argument: contraction. The parameter 6 was chosen such that 
the scaling parameter a{t) belongs to a fixed compact subset of (0, oo) and therefore the 
imaginary eigenvalue ia fulfills 

G [01,02] C (0,00), (2.107) 

for all the admissible paths that we consider. Fix, then, a constant p £ (0, ai). 

For any two bounded solutions of the linearized equation (|2.33p . {Rj,7rj) = T(i?°,7r^), 
j = 1,2, located in X, such that 

\\{R^j,7r^)\\x<d, (2.108) 
we prove that T acts as a contraction in the following space Y: 

Y = {iR,7T) \ \\e~'PRit)\\^^^r/,^^,^^;,,e + \\e-'Pn{t)\\Li < 00}. (2.109) 

This is only a seminorm, but defines a metric space for fixed initial data. 
Furthermore, for fixed Rq we prove that the unique value of h, H^Rq, R^,7r^), for which 
the solution with initial data (j2.56p is bounded, satisfies 

|Mi?o,i??,vr?) - /i(i?o,i?^,vrO)| < C<5||(i??, vr?) - (i?^, 7r0)||y. (2.110) 

Observe that this enough to complete the proof. For initial data given by ()2.56p and 
fixed initial data Rq and ttq, consider a sequence 

{Rn,nn) = T((i2„_i,7r„_i)), ||(i?„, 7r„)|U < 5 (2.111) 

which converges in the Y sense to (-R, vr); the parameters /i„ = h{Ro, R^,Tr^) form a 
Cauchy sequence as well. Then the pair (R, vr) is a fixed point of T and, by virtue of 
Lemma 12. 2| a solution to the nonlinear equation (locally in time in a weak sense and 
therefore globally as well) with the specified initial data and, furthermore, 

||(i?,vr)|U <limsup||(i?„,7r,)||x <<5. (2.112) 

This follows first on any finite time interval [0, T] and then in the limit on [0, cxd). 
We seek to prove that for any sufficiently small choice of 5 

||(i?i,7ri) - (i?2,vr2)||y < 1/2|| (i??, vr?) - (/?°, ^2°)l|y • (2.113) 

Letting Ri and R2 start from distinct initial data proves useful, leading to the formu- 
lation of the following perturbation lemma, which we employ repeatedly: 
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Lemma 2.5. Consider two solutions, (i?i,7ri) and {R2^'^2), of two distinct linearized 
equations in the form h2.2>2>\) : 

idtRj + n,o{t)Rj = Fj, Fj = -iL^oRj + iV(i?o, w^o) - iV,o(i?o, W^o) 

f, = 4a'^\\W^o\\f{{R,,id^EfiW^o))ir'^) (2.114) 

f, = 2\\W^o\\^\{R,,id^Ef{W^o))ir^) -i{N{RlW^o),E^^ {vk,Dk}. 

for j = 1,2, with initial data 

Rj{0) = Roj + hjF-^iWiTro)), TTj{0) = 7r°(0) = ttq given. (2.115) 

Assume in addition that || (ii^, 7r^)||x < (5 and 

{Rj,7T,) = T((i?o, ttO)), hj = h{Roj,R°j, vrO), (2.116) 

meaning that hj hold the unique values that make {Rj,TTj) the bounded solutions of ( 12.114|] . 
Then, assuming 6 > is sufficiently small, 

||(i?i,7ri) - {R2,7r2)\\Y < C^IlK,^?) - {R^,4)\\y + C\\Roi - i?02||^i/2, 
- /12I < C5(||(i??,7rO) - (fl°,7rO)||y + ||i2oi - ^02^1/2). 

In this formulation, 6 and various constants denoted by C may depend on the value of 
p in the definition of Y, (j2.109p . which in turn depends on the scahng parameter a. 

Proof. Let R = Ri — R2, = tti — tt2- Rj, j = 1,2, satisfy the equations 

idtRj + n^o{t)Rj = Fj, (2.118) 

with initial data 

Rj{0) = Roj + hjF+{W{Tro)). (2.119) 

Note that, by the preceding stability result. Proposition 12.31 IK-^j; < ^ for j = 1,2. 

Furthermore, Rj satisfy the orthogonality conditions 

{R,it),EjiW^om=0 (2.120) 

at time t = and thus, due to equation (j2.114p . at every time t. 

Subtracting the linearized equations from one another, we obtain a similar one for the 
difference R = Ri — R2: 

idtR + n^o{t)R = F, F = Fi-F2- {n^o{t)-n^o{t))R2. (2.121) 

We choose the Hamiltonian Tl^o (t) (the choice of one or two is arbitrary) and apply the 
same isometry U{t) as described by (12. 65112. 681) : 

U{t) = e/o + ^((«?)'(^) - K{s)\')as) ds _ (2.122) 

Let 

Z{t) = U(t)R{t). (2.123) 
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We reintroduce notations similar to ()2.74p . namely 



m 




Po{t) 


= Po(T^(vr?(i))), 


Pc{t) 


= Pe(T^(vr?(i))), 


Pimit) 






= p+(H^K(t))), 


P-{t) 


= P_(m<(i))), 




= F+{w{4m, 


F-{t) 








a{t) 





(2.124) 



Rewritten for Z, the equation becomes 

idtZ + n{t)Z = U{t)F. (2.125) 

Then we split Z into three parts, according to the Hamiltonian's spectrum: 

Z = P,{t)Z + PQ{t)Z + Pi^{t)Z. (2.126) 

We bound each component of Z in the same manner as in the previous section. The 
main difference is related to Po{t), which need no longer be zero because the orthogonality 
condition does not hold. 

For the right-hand side F = Fi — F2 — {TC^o{t) — Ti.^o{t))R2 we estimate the difference 
term by term, arriving at the following estimate: 

l|i^lL*Pi2^iA6/5 < C6{1 + p-2)(||(i?,vr)||y + \\{RUi) - (^2,^2)l|y)- (2.127) 

Here we have taken advantage of the exponential weight, as integrating in time preserves 
the space e^^L\, at the cost of a factor of : 

\\!lf{s)ds\l,^^,<Cp-^\\f\U,j^v. (2.128) 

In particular, one has in any Schwartz seminorm 5^, in the space variables, that 

\\W^o{t) - W^.{t)U^s. < C(l + p-'){\\{Rl4) - {RIA)\\y (2.129) 

The initial data are given by 

Z(0) = {hi - h2)F+{W{^o)) + Roi - R02 (2.130) 

and satisfy the bound 

11^(0)11^1/2 < C(|/ii - /12I + ||i?oi - ^0211^1/2)- (2.131) 

It is straightforward to bound the continuous spectrum projection, Pc{t)Z, which sat- 
isfies the equation 

idtiPcit)Z) + n{t)P,{t)Z = Pc{t)U{t)F + i{dtPc{t))Z. (2.132) 

Starting from the usual Strichartz estimates of Theorem 13.91 and integrating against the 
exponential weight e*'', we obtain Strichartz estimates with this weight: 

ll^c(t)^||^,p(^^^^l/2nL2vt.l/2.6/5) 

<C(\\PMZm\m/^ + ll^lle.PL2^iA6 + \\dtPc{t))Z\\^,^^,^,j,,,) (2.133) 
<C||i?oi -i?02|lHi/2 +C'^(1 +P~')(ll(^?,0 - (^2,^2)l|y + ll(^,vr)||y). 

There is no contribution due to hi — /i2, since Pc(0)-^(0) does not depend on hi — h2- 
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R satisfies no orthogonality condition, but Ri and R2 do. Taking the difference of these 
two relations leads to 

(i?i(t)-i?2(t),H^(W^.o(t))) = (i?2(t),H;(W^,o(t))-H;(M^,o(t))). (2.134) 

Applying the isometry U, it follows that 

{Zit),Efit)) = {R,it)-R2{t),Ef{W^oit))) (2.135) 

and therefore 

||Po(t)^(t)||,,,(^o.^i/2n^2,^iA6/5) < C5(l + /;-2)||(ii?,7r?) - {R^2,4)\\y. (2.136) 

Concerning the imaginary component, Lemma 12.41 applies again, since Z is globally 
bounded forward in time, as seen from 

\\Pimit)Z{t)\\^,/, < \\Rimm/2 + \\R2mm/2 < C < 00. (2.137) 

Write, following Section [2T2l 

P^mit)Z{t) = b+it)F+{t) + b-{t)F-{t). (2.138) 

Then, as in the proof of the stability result, we infer 

dtb^ = ±a{t)b^ - {UF,asF^)- 

- 3d°(a°)-^(Z,m3FT) + (a°)-=^(Z,m3((i^FT)^°). 

Thus b^ and b~ satisfy the equation 

where 

N± = -{UF,asF^)-3a'i{a'i)-^{Z,iasF^) + (a?)-3(Z, m3(d,FT)0. (2.141) 
Then, 

\mL.,L^<C6\\{R,7r)\\Y. (2.142) 

As a consequence of the lemma, we obtain that Pim{t)Z is in the same space as A'^-i-, as a 
function of time. Indeed, here we emphasize that since p < ai < cr(t), convolution with 
g-«il*l preserves the space e^^L^. One gets that 

("OO 

\hi - h2\ = \b+{0)\ < C / e-*"i|A^+(t)|dt 

Jo (2.143) 

< C<5(||(Z?,vr?) - (ZO,vrO)||y + 

Also through Lemma 12.41 it follows that 

mmZ\\Y < C<5(||(Z?,vr?) - (Z0,7r0)||y + ||(i?, 7r)||y). (2.144) 

Finally, subtracting the modulation equations (j2.114p . the difference between the paths, 
vr = TTi — 7r2, also fulfills 

Me^.Ll < C5(l + p-2)(||(ZO,vr?) - {zI,4)\\y + ||(Z,vr)||y). (2.145) 

Up to this point, there was no actual need for exponential weights and we could have used 
polynomial weights in their stead. However, (|2.145p makes exponential weights necessary. 
The modulation equations involve terms of the form / on the left-hand side and Jq f{s) ds 
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on the right-hand side. The only way to contain both within the same space involves 
exponential weights. 

Putting (|2.133p . (|2.136|) . (|2.144l) . and (I2.145|l together we retrieve estimate (|2.117|1 . 

□ 

Remark 2.6. A difference estimate also exists in the case of different starting solitons, 
7r°(0) = 7rj(0) = TTQj, j = 1,2. In order to perform a comparison, we first make the starting 
solitons coincide by means of a symmetry transformation and then compare the solutions 
in a unified setting. 
More precisely, let 

Q{t)f{x, t) = e^(''-^-*l''l'+^)a/(ax - 2tv - D, ah) (2.146) 

be the unique symmetry transformation such that g(0) takes w{ttoi) to w{'itq2). We define 
accordingly the action of g on the modulation parameters: 

QTT = TT Qw{7r) = w{tt). (2.147) 

Then, using the notation q{R,tt) = (g(t)i?(t), 0(t)7r(t)), it immediately follows from the 
previous lemma that 



\e{Ri,7ri) - (i?2,^2)||y < C5||0(i??,7r?) - (i?o, 7r0)||y + C||s(0)i?oi - i?02||^ 

\hl - h2\ < CSiMR^TT^,) - (i?^,^2°)l|y + M0)R01 - R02\\hi/2] 



(2.148) 



2.5. Analyticity of the invariant manifold. Consider a solution of the nonlinear equa- 
tion (jO.ip stemming from the contraction argument presented. For a soliton Wq = W{'Kq), 
it has the form 

i? = Wo + i?o + KR^, Wo)F+(Wo). (2.149) 

Here F+(VFo) is the eigenvector corresponding to the upper half-plane eigenvalue of TLiyVo) 
(see Section [2T2|) ■ Rq belongs to the codimension-nine vector space 

A/'o(VFo) = (Pc(Wo) + P-{Wo))H^'\ (2.150) 

and /i(i?05^o) is the unique value determined by the contraction argument that leads to 
an asymptotically stable solution to (10. Ih . for these initial data. 

We recall again that R and other vector-valued functions we are considering have the 

form ( ) and that, as a consequence, the dot product is real-valued. 



r 

At this point we give the following formal definition: 

Definition 2.2. Let M be the eight- dimensional soliton manifold and 

Mo{W) = {R^ G {Pc{W) + P.{W))H^/^ I \\Ro\\hi/2 < So} 
Af{W) = {W + Ro + h{Ro, W)F+{W) \ Rq G Mo{W)} 

Afo = {{Ro, W)\Ro€ Afo{W), W e M} (2-151) 

Af= \J Af{W). 
5q will be chosen independently of W. 
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The fiber bundle A/q is trivial over the soliton manifold. Indeed, for each soliton W 
there exists a unique symmetry transformation qw that takes Wq (a fixed soliton) into 
W. Then 

{Ro,W) ^ QwiWo + Ro) (2.152) 
is an isomorphism between a tubular neighborhood of the base in the (trivial) product 
bundle {Pc{Wo) + P-{Wo))H^^'^ x M (where M is the soliton manifold) and A/q. This 
endows A/q with a real analytic manifold structure. 
M is the image of A/q under the map 

J'iRo, W) = W + Ro + h{Ro, W)F+{W). (2.153) 

Following the contraction argument from beginning to end and giving appropriate values 
to 6, commensurate with the size of the initial data Rq, we summarize the conclusion as 
follows (this is the gist of our main theorem): 

Proposition 2.7. There exists 5o > such that for each soliton Wq = W{7rQ) there is a 
map h{-, Wo) : AAo(Wo) ^ ^ such that 

(1) h is locally Lipschitz continuous in both variables, 

(2) |/i(i2o) W^o)! ^ C'ao 11-^0 11^1/2 (Q^o is the scaling parameter ofWo), 

and 

HRo, Wo) = Wo + Ro + h{Ro, Wo)F+{Wo) (2.154) 
gives rise to an asymptotically stable solution ^' to i0.1\) with ^'(0) = J^{Ro, Wq) such that 

^{t) = W^{t) + R{t). (2.155) 

Here M^- is a moving soliton with W,r(0) = Wq, governed in accordance to il.2A]) by a 
path vr such that 

Ml < Cao\\Ro\\%^f2- (2.156) 
R is in the Strichartz space, with initial data R{0) = Rq + h{Ro, Wo)F^ (Wo) , and 

Adopting a different point of view, for sufficiently small So \\Ro\\^i/2 is comparable to 
the distance from ^'(0) to the soliton manifold, 

rmnll^'(O) - (2.158) 

Indeed, one inequality is obvious and the other follows by Lemma 12.121 see below. 

Applying the perturbation Lemma 12.51 to the solution of the nonlinear equation leads 
to the following: 

Proposition 2.8. The solution depends continuously on the initial data: given solutions 

^j{t) = W^^{t)+R,{t), (2.159) 
j = 1, 2, with initial data in the same fiber, \I'j(0) G M{Wo), meaning ^'j(O) = J^{Roj, Wq), 



^1 - ^2 ll,.p(ic„ ^1/2^^2^1/2,6/5) < C||(i?i,7ri) - (i?2,vr2)||y ^^^^ 



one has 

11^1 - ^2 ||^,p(^c„ ^1/2^^2^1/2.6/5) 

< C||iioi — -R02II//1/2 
and 

\h{Roi,Wo) - h{Ro2,Wo)\ < Cao ( II i?oi 11^1/2 + ||i?o2bi/2)||i?oi - i?02 11^1/2- (2.161) 
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The continuous correspondence between different fibers is given by symmetry transfor- 
mations. 

Proof. Firstly, applying the comparison Lemma l2 . 51 yields 

||(i?i,7ri) - (i?2,vr2)||y < C||i?oi - i?02bi/2. (2.162) 

The comparison Lemma 12.51 also implies 

\hiRoi,Wi)-h{Ro2,W2)\ < CS\\Roi- Ro2\\hi/2. (2.163) 

Taking 6 proportional to ||-Roi |li/i/2 + II-R02 ll^fi/2 leads to the second conclusion. □ 

In particular, this shows that the map given by (j2.153p is locally Lipschitz continuous. 
We explore its properties further, beginning with a definition and a preliminary lemma. 

Definition 2.3. Given two Banach spaces A and B, a map f : A ^ B is analytic if it 
admits a Taylor series expansion: 

f{a) = fo + fi{a) + /2(a,a) + . . . , (2.164) 

such that for each n fn is n-linear and there exist constants Ci, C2 such that ||/n|U®"->_B ^ 
CiC^. 

A similar definition can be given for differentiable, C" class, and smooth maps. Fur- 
thermore, the definition naturally extends to manifolds. 

Lemma 2.9. The map F : TVq x M ^ H^l'^ , 

F{W,R,h) = W + R + hF+{W), (2.165) 

is locally a real analytic diffeomorphism in the neighborhood of each point (VF, i?, 0). 

Proof. Let Wq = W{'Kq) and consider the differential of -F, given by the linear map 
dF \(WM,y- X M)(W^o) X M ^ H^^ 

dF \iWo,RoM) {^T^.^R, Sh) = {d^Wo)67T + 6R+ {6h)F+{Wo) + H{d^F+ {W^))8^. 

(2.166) 

That dF is bijective at points where /iq = follows from the following identity: 
= Pq{Wo)^ + Pc{Wq)-^ + P-{Wq)-^ + P+(Wo)^ 

= {^.^f{Wo))dfWo + {Pc{Wo) + P-{Wo))^+ ,2 167) 

/e{a,r,t;fc,Dfe} 

+ a^''{^,icj^F~{Wo))F+{Wo). 
This leads to an explicitly constructed inverse for the linear map: if 

^ = dF \(^v,Mo) {S7T,6R,6h), (2.168) 

then 

67r = ((VI/, EM)), ^t{Wo)), H,, (I^o)>, ^d, (Wo))), 

5R = (PciWo) + P-iWo))^, (2.169) 

6h = a^^{^,ia3F-{Wo)). 

The local invertibility of the nonlinear map F follows by the inverse function theorem. 
Smoothness follows by inspection of the explicit forms of dfWo and dfF^{Wo). 
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Next, we consider the analyticity of which is closely tied to that of the soliton 
W{'7t), of its derivatives, and of the eigenvectors F~^{W{Tr)), considered as functions of the 
parameters vr. 

The analyticity of the soliton was shown by Li-Bona [LiBoj . but we need a stronger 
statement than the one proved in their paper. 

Lemma 2.10. Given an exponentially decaying solution (j) of equation /i0.3\) 

- A(j){-, a) + a'^(j){-,a) = a), 

the soliton 

w{7r){x) = e^(^-^+^V(2; - D, a) (2.170) 
is a real analytic function of it in any Schwartz class seminorm S. Furthermore, 

||5^a*5t9g'=z/;||5<CiCfl \{ (3fl. (2.171) 

We distinguish between two kinds of parameters in the definition of vr (j2.170p . In the 
parameters a and v^, the proof makes it clear that the domain of analyticity is not the 
whole complex plane, but only a strip along the real axis. In the parameters F and Df., 
on the other hand, the soliton extends to an analytic function of exponential type in the 
whole complex plane. 

Proof. To begin with, we show that w is an analytic e~^^~^^l^lL^-valued map. Since the 
derivatives of an analytic map are also analytic and dc^'^ij^) = —dx^,w{7r), by iterating it 
follows that arbitrarily many derivatives of w^tt) are analytic maps into this space. 

The real analyticity of w is equivalent to its complex analyticity and the joint complex 
analyticity is equivalent to separate analyticity in each variable. Given that 4>{-,a) decays 
exponentially at the rate e~^" by the Agmon bound, we can extend w to an analytic 

function for all T and on the strip {v \ \ Imf | < a^}. This proves analyticity in regard to 
r and V. 

Analyticity in the other two variables, D and a, requires that we show that (/>(-, 1) is an 
analytic function in the spatial variables, in the sense that there exist constants Ci and 
C2 such that for every multiindex (3 = {(3^-^ , (3^^ , /3t,3 ) 



12 <Ci|/5|!C^'^'. (2.172) 
4> is exponentially decaying, |i;^(a;)| < Ce"^^"*^-"^', and satisfies the equation 

-A0 + (/> = g((/>), (2.173) 
where g{z) = is an analytic function. Differentiating (3 times, we obtain 
{-A + l){d'^(P) = d>^igocP) 

m 

= E E d^^cl>-...-d^^<P-ig^^Kcl>). 

3=1 Pi+...+l3j=p 



(2.174) 



In the proof we make the more general assumption that g is analytic and its derivatives 
grow subexponentially, meaning that there for any C4 > there exists C3 such that 

\9^'\^)\<C,Cl. (2.175) 
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It follows that for any C4 > there exists C3 such that, for fixed / G i^^, 

II^Wo/ll^, <C73Ci. (2.176) 

This property characterizes not only polynomials such as the case of interest, g{z) = , 
but also functions of subexponential growth, such as g{z) = cosh(-y/z). 

To avoid complications, we present the argument in the algebra H'^. We prove by 
induction on \j3\ that, for properly chosen Ci and C2, 

Wd'Hu^ < ^fl- (2.177) 

Indeed, assuming that this induction hypothesis holds for all indices up to |/3|, we obtain 
that 

\^,<c\\{-A+i){dfmH^ 

<CCd^\y^ V gjgf (2.178) 

<^^3C2 2^ 2^ (|/3i|2 + i).....(|/3J2 + i)- 



Observe that 

n 

^ (fc2 + l)((n - fc)2 + 1) " n2 + l 



(-^2 4_ 1 VCr, - _n ^1 < „2 _n • (2.179) 

fc=0 



By induction, we obtain that 



1 6^-^ 

> < -. 2.180 

^ nf + 1) • . . . • (n2 + 1) n2 + 1 



ni-|-...-|-nj=?i ^ J- ' ^ J 

Therefore 

By making C4 sufficiently small and C2 sufficiently large, we obtain that the sum is 
uniformly bounded, regardless of |/?|, and 

This proves that the induction assumption also holds for derivatives of order + 1. 

The preceding proof shows that 0(- — D) is an analytic -ff^-valued map. However, the 
argument works for any algebra A containing and with the property that 

||d/|U<C||(-A + l)/|U. (2.183) 

This means that the algebra A has to be invariant under convolution with the kernels 

^ + (2.184) 

Note that the algebra A = e^'^^^'^)\^\lJ^ fulfills both requirements. Thus, 0(- — D) is an 
^-valued analytic map. This ensures the joint analyticity of w{ti^ given by ()2.170p in the 
variables -D, f, and F around the point ^ = ii;(l, 0, 0, 0). By symmetry transformations, 
this implies analyticity around any other point. 
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Concerning a, observe that the derivative with respect to a is given by the generator 
of dilations, which is a combination of multiphcation by x and taking the gradient V: 

da(l> = a-^(f) + xV(p. (2.185) 

Therefore, analyticity in a follows from that with respect to to D and v. 

If W is analytic, so are its derivatives up to any finite order, enabling us to conclude 
that W is analytic in any Schwartz class seminorm. □ 

In a similar manner, the analyticity of (W {tt)) as a function of vr reduces to that of 

//+\ 

a fixed eigenfunction F+. = ( j+ 1 (see Section [2T2]) satisfies the equation 

(A-l + 2,/.2)/+ + </.27+ = ia/+. (2.186) 

Knowing that (p is analytic, proving the analyticity of proceeds in exactly the same 
manner as above. 

Clearly, if W{tt) is analytic then so are its derivatives up to any finite order. The 
same goes for F^{W{7r)) and its derivatives. This immediately implies that F is real 
analytic. □ 

Lemma 12.91 has the following immediate consequence: 

Proposition 2.11. given by 12.153)) is locally one-to-one and its inverse (defined on 
its range) is locally Lipschitz. 

Proof. The local invertibility of J- follows immediately from the previous lemma. Indeed, 
one has 

T{R,W) = F{W,R,h{R,W)). (2.187) 

For a sufficiently small Sq, h{R, W) is close to zero and the previous lemma applies. In order 
to establish the Lipschitz property for the inverse, we can simply ignore the parameter 
h. □ 

Another consequence is that, if a function is sufficiently close to the manifold Ai of 
solitons, we can project it on the manifold as follows. 

Lemma 2.12. For every soliton W there exists 5 > such that whenever \\^ — W\\jji/2 < S 
there exists Wi such that Pq{Wi){^ — Wi) =0 and 

11^ - Wi\\^i/2 < C\\^ - (2.188) 
Furthermore, Wi depends Lipschitz continuously on ^. 

Again, following the use of symmetry transformations, S can be chosen without regard 
for W. 

Proof. If 6 is sufficiently small, ^ = F{Wi,R,h) for some {Wi,R,h) close to (1^,0,0) by 
Lemma 12.101 Since the inverse of F is bounded, 

\\W - Will^i/a + ||i?||^i/2 + \h\ < C\\^ - W||^i/2 < C6. (2.189) 

Since R G Mo{Wi), by definition Po(W^i)(^ - Wi) = 0. 

To a first order W — Wi lies in a direction tangent to the soliton manifold M , meaning 
within the range of Po{Wi), so by means of a Taylor expansion we obtain 

||(/ - Po{Wi)){W - PFi)||^v2 < C\\W - WiWl,,,. (2.190) 
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Thus 

< C||^ - W\\^,/, + - Po{W{)){W - Wi)\\^,/2 (2.191) 
<C(||^'-PF||^,/2 + 11^-^^111^,1/2)- 

On the other hand, 



(2.192) 



11^^ - Will^i/2 < (11^ - W^lli,i/2 + 11^- - WiW^.nf 

< C\\^-W\\^,^2+C6\\^-Wi\\^,^2. 
For 6 sufficiently small, the conclusion follows. □ 



This lemma would be superfluous if Pq were an orthogonal projection and the constant 
could be taken to be one then. Also note that in this generality the conclusion still holds 
for VF^/-^'^, with the same proof. 

Definition 2.4. By small asymptotically stable solution we mean one that can be written 
as ^{t) = M^7r(t) + R{t) where WV(t) is a moving soliton governed by the parameter path 
TT as in i fl.24)) and 

||(i?,7r)||x = Pllioo ^1/2^^2^1/1/2,6 + ll^lli < 5o. (2.193) 



X is the space that appears in the contraction argument, see (I2.59j) . 
We can rewrite any small asymptotically stable solution ^ as Wj^{t) + R{t) such that 
the orthogonality condition is satisfied: 

Po{Wmt)))R{t) = 0. (2.194) 

Following Lemma |2.12| R{t) is still small in the space 

Furthermore, VKs-(t) depends Lipschitz continuously on ^. Writing the modulation 
equations explicitly as in ()2.2ip . it follows that ||vf||/^i is small too. 

Thus, it makes no difference whether we assume the orthogonality condition initially, 
as part of Definition 12. 4| since we can produce it in this manner. We arrive at the same 
definition with or without the orthogonality condition, though it may be for different 
values of 60 . 

Clearly, every solution with initial data on the manifold M is small and asymptotically 
stable. A partial converse is also true. 

Proposition 2.13. If^{0) is the initial value of a small asymptotically stable solution ^ 
in the sense of Definition \2M to equation i0.1\) . then ^'(0) € M. 

Proof. Write ^ = WT^^t) + R{t), satisfying the orthogonality condition 

Po{W^{t))R{t) = 0. (2.195) 

By Lemma [2. 91 there exist Rq £ A/'o(VF7r(0)) and h such that 

R{0) = Ro + hoF-^iW^O)). (2.196) 

Then, we note that both the initial data 

*(0) = W^{0) +Ro + hoF+{W^{0)) (2.197) 
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and the initial data 

^{0)=W^{0) + Ro + h{Ro,W^{0))F+{W^{0)) (2.198) 
give rise to small asymptotically stable solutions, in the form 

^{t) = W^{t) + R{t), ^{t) = W^it) + Rit). (2.199) 
The perturbation Lemma 12.51 then applies, implying that 

||(i?,7r) - {R,n)\\Y < C6\\{R,7r) - (i?,*)||y. (2.200) 

Otherwise put, we obtain that {R,it) = (i?, vf). Applying the lemma once more, it follows 
that the h values coincide as well, that is ho = /i(i?Oi W^7r(0)). Therefore 

^i>{0) = ^{0)=W^{0) + Ro + h{Ro,W^{0))F+{W^{0)) (2.201) 

and thus ^'(0) belongs to M. □ 

Corollary 2.14. If ^ is a solution to Ii0.1\) whose initial data ^'(0) belongs to N , then 
^{t) also belongs to N for all positive t and for sufficiently small negative t. 

Proof. Clearly, both for positive t and for sufficiently small negative t ^'(t) exists (due to 
the local existence theory, for negative t) and gives rise to a small asymptotically stable 
solution. Then, the previous proposition shows that ^{t) must still be on the manifold. □ 

To recapitulate, we have investigated the properties of five maps that describe the 



solution of the nonlinear problem (jO.ip , namely 

/i(i?o,Wo) : AAo ^M, (2.202) 

T{Ro, Wo) :Mo^J\f, T{Ro, Wo) = F{Wo, Ro, h{Ro, Wo)), (2.203) 

and the solution itself, 

^ = ^(i?o, Wo) : AAo ^ n LfW^/^'^/\ "^{t) = W^{t) + R{t), (2.204) 

TT = Tr{Ro, Wo) : Afo ^ Wl'\ (2.205) 

R = R{Ro, Wo) : Mo L'^H^^^ n L^Vi/^.e/s^ ^2.206) 



It turns out that all five are real analytic and one can also replace A/q by — change the 
variable to ^'(0), that is. 

Proposition 2.15. For fixed Wo there exists p > such that the maps h, T considered 
as a map into e^^Lf , (i?, vr) regarded as a map into Y (see < [2.109|) ). and seen as 
map into e^f{L^Hx^'^r\L^Wx^'^'^), are real analytic in the variable Ro on the fibre Mo{Wo) 
or, equally, in the variable ^'(0) on the fibre M (Wq) ■ 

Letting Wo vary, note that the dependence of these maps on Wo (across the different 
fibers of A/q) is given by symmetry transformations, which are also analytic in a properly 
considered setting. 

Also note that the analytic dependence of the solution on initial data and the analyticity 
of the manifold are closely tied to that of the nonlinearity in (jO.ip . If the nonlinearity 
were of class C'^, we would expect that the manifold and the dependence on initial data 
would also be of class C^. 
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Proof. Firstly note that, since F is analytic, it suffices to show that h and vr) are 
analytic. 

The proof consists in the following steps: firstly we exhibit the first-order differential of 
h and (i?, vr) and show that it satisfies the definition of differentiability, meaning that the 
remainder in the Taylor expansion is quadratic in size. By recursion, we then introduce 
the n-th order differentials for arbitrary n. 

Finally, we show that differentials and remainders in the Taylor expansion grow in norm 
at most exponentially, enabling us to conclude that h and {R, vr) equal the sum of their 
own Taylor series within a sufficiently small convergence radius. 

To begin with, for given initial data Wq = W{ttq) and Rq, denote the corresponding 
solution to (jO.ip by ^' = Wj^ + R. R and vr satisfy the equation system, in the form (j2.33p . 

idtR + n^{t)R = F{Tr,R), 

(2.207) 

f = Ff{7T,R),fe{a,r,Vk,Dk}, 

where 

, ! , , , , 2.208) 

F{Ti,R) = -iL^R + N{R,W^)-N^{R,W^), ^ ' 

and we introduced the notation 

' Aa\\W^\\'^^{{R, {d^Ef{W^))ii) 

-i{N{R,W^),~f{W^))), /G{Q,r} 
2\\W^\\-^^{{R, KH/(Ty^))7^) 

-i{N{R,W^),~j{W^))), fe{vk,Dk}. 
The initial data are given by 

i?(0) = R^ + h{Ro, Wo)F+{Wo), 7r(0) = ^o- (2.210) 
This setup makes the orthogonality condition 

{Rit),Ef{W^m = (2.211) 

valid for all times. Indeed, it holds initially at time t = and equation (I2.207P ensures 
that it still holds at any other time. We aim to find an infinite expansion of the form 

i? = i?o + i2i + i?2^..., vr = 7r° + 7r^ + 7r2 + ..., h{Ro,Wo) = h° + + + . . . , (2.212) 

where i?" = n\d"'R, vr" = uIcPtt, and /i" = nld^h are the n-th order terms in the power 
expansions of R, vr, and h, given by n-linear expressions in i?o, such that 



Ff{7r,R) 



(2.209) 



\R- R'^ - R^ - ^"lletp(ico^i/2nL2iyi/2'6) ^ <^™||i?o||^t/2, 



_ ^0 _ ^1 _ ^ ^ _ _ ^n||^^^^^^ < Cn\\Ro\\f,%, (2.213) 

Showing, in addition, that C„ < CiCg for some constants Ci, C2 and every n guarantees 
the analyticity of R, vr, and h. 

For technical reasons, it is easier to construct the differentials and prove analyticity at 
zero. Analyticity will then follow at every point within the radius of convergence. 

The constant terms in the expansion are 

RO{t) = 0, 7r°(t) = vro = (1,0,0,0), /i° = 0. (2.214) 
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The natural guess is that first-order differentials and vr^ satisfy the following lin- 
earized version of (j2.207p : 

idtR^ +n^o{t)R^ = (a^F(^°,i?°))7ri + {dRF{TT^,R^))R^ - (a^Ko)(t)vr^i?° 

(2.215) 

/I = {d^Ff{n^,R^))7T^ + {dRFj{7T^,R^))R\ f G {a,T,Vk,Dk}, 
with initial data 

i?i(0) = Ro + h^F+{Wo), iT^{0) = 0. (2.216) 

It is important to note that, because i?*' = 0, vr'^ = 0, the following terms cancel: 

id^V^o){t)7T^R^ = 0, d^Ffin"^, i?°) = 0, dRFfin'', i?°) = 0, 

5^F(^o, i?0) = 0, dnFin^, i?") = 0. ^^'^^ 

This improvement only holds when taking the differential at zero and is due to the fact 
that the nonlinearity is of order higher than two. 
It follows that TT^ = and 

idtR^ + n^o it)R^ = 0. (2.218) 

The orthogonality condition 

{R,Ef{W^)) = (2.219) 

becomes 

{R\Ef{W^o)) + {R',{d^Ef{W^o))7T') = 0; (2.220) 
consequently it still holds for R^: 

{R\EfiW^o))=0. (2.221) 

By means of Strichartz estimates we obtain, in the same manner as in the proof of stability, 
that R^ is bounded for a unique value of and then it satisfies 



\^^h^Hl/^nL^,w^^^'^ - (2.222) 



Note that rcx) ui/2„r2Ti/i/2,6 cau be made arbitrarily small, but ||7T^||ri = 0. 

To show that R^ thus constructed and vr^ = are indeed the first-order differentials of 
R and vr, consider 

= R^ + R\^^ = TT^ + tt\ (2.223) 
If F{tt,R), Vj^, and Ff(TT,R) were differentiable, it should be the case that 

F{^\S') = F(i?o, vr") + (a^F(^o, R^))7t' + {dnF{7r^, R^))R' + oi(i?\ tt^), 

Vj^i = V^o + d^V^oTT^ + 02 (^'), (2.224) 
Ff{^\S') = F^(vrO,i?0) + (a^F^(vrO,i?0))^i + {dRFfi7r'',R^))R^ + 03(i?\^i). 

In fact, even more is true: all three quantities are analytic, as one sees by examining their 
explicit forms. 

The main issue is that the soliton W7r(t) depends not only on the values of vr and vr at 
time t, but also on the integral Jg vr(s) ds. Thus, even though W^n- depends analytically on 
vr, each derivative produces a factor of t. 
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In this setting, the error terms in (j2.224p are quadratic: 

ll^2(7r^)||(t^^oo(iyi/2,6/5-.niyi/2,6/5+.) < C 1 1 (i?^ TT^ ) 1 1 ^ (2.225) 
\\os{R\7r')\\^,^L.<C\\{R\7r')\\l. 

Therefore 



fj,r=Ff{^\S')+0}iR\7r'), f G{a,T,Vk,Dk}, 



(2.226) 



where O^, are error terms bounded by 

||02(i?i,vri)||^^^^,^V..a/5 + \\0}{R\7r%,^L^ < C\\{R\n')\\j,. (2.227) 

Following (j2.22ip . fulfills the orthogonality relation 

{S\t),Ef{Wj,im = 0. (2.228) 

Comparing the equation system (j2.226p satisfied by and with the one satisfied by 
R and vr, (|2.207p . we obtain that 

||(i?,vr) - {S\j:')\\y < C\\{R\7r')fj, < C\\Ro\\l,,,. (2.229) 

The proof exactly follows that of Lemma 12.51 

We repeat this procedure for the higher-order terms in the expansion. Let n > 1 and 
consider a variation it = ir^ + Sir of tt" such that (57r(0) = 0. Recall that solitons depend 
analytically on parameters: for any Schwartz seminorm S, 

||d^VF(7r°)(5^))||s < CiC^n\\\d^'l. (2.230) 

However, by ()2.17ip only derivatives in a and produce a factorial contribution: 

/e{a,-Ufc} 

Thus we obtain for W^o given by (fOil) and ([230]l 

w^oit) = w(^a^{t)y{t) + {{a°{s)f - \v^{s)\'^) ds,v^ (t), D{t) + 2 v°{s)ds 

— and likewise for djW^o, F~^{W^o), and all other quantities that depend on the moving 
soliton — the following explicit expression of analyticity: 

||a^a*a^-9^^^VF.o(t)(<5vr))||5<CiCfl n ^f'- n W^'II^^IlP- (2-231) 

As noted before, this is an improvement over mere analyticity (and one that is necessary 
in the sequel), in that, while in half the variables W is only analytic on a strip, in the 
other half W extends to an analytic function of exponential type on the whole complex 
plane. 
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At our discretion, we pick parameters a and ai such that a < ai < a{a{t)), for any time 
t, where a{t) is the scahng component of 7r(i) and vr is any path that appears in this proof 
(their scahng parameters are uniformly bounded away from zero). Define the weights 



^„(t) = yi^<e(«*>. (2.232) 
j=o 

Then imphes 

W^o = W^o + Wlo + ... (2.233) 
such that in any Schwartz seminorm Sx, in the space variables only, 

ll^;o(5vr)|U„(,)ic«s^ < CCl^Wdirr^. (2.234) 
This immediately yields a power series expansion for V^o: 

= V^, +Vl, + ..., (2.235) 

with V^n n-linear for each n and 

\\y:o{87:)\\A^it)LTs. < cicnsnr,. (2.236) 

The weights An{t) given by (j2.232p have the property that 

n 

Mt)An-j{t) < C7M„(t). (2.237) 

j=0 

Consequently, consider the spaces d^^An{t)Ll. In this setting, for 

mi + . . . + nin = m, (2.238) 

and given n path variations Stti to ^vr^ such that Sttj = 0, j = l,n, one has that 

W {67:1, ... ,6nn)\\A^+„(t)L^S^ < CiC2"||(57^i|U^^(t)ii • ... • \\SnJA^^(t)Lj- (2.239) 

Considering the explicit form of F{R,Tr) and Ff{R,iT), note again that 

i^(i?°, 7r°) = dF(i?o,7r°) = 0, F/(i?°,7r°) = dF/(i?o,7r°) = 0. (2.240) 

This is due to linearizing around zero {R^ = 0, tt'^ = 0), starting from a nonlinearity of 
order higher than two. 

In fact, one needs to take two derivatives in i? or vr before arriving at a nonzero derivative 
of F OT Ff. We obtain the power series expansion 

F = F° + + . . . , Ff = Ff + F) + (2.241) 
where F" and Fj? are n-linear and, for X defined by (j2.59p . 

II^"IIa .,.^2^1/2,6/5 <CiC2"||(5i?i,<57ri)|U-...-||(5i?n,<5^„)|U, 
IU„_,(t)Li < CiC2"||(<5i?i,57ri)|U • . . . • ||(5ii„,<57r„)|U. 

The gain from n to n — 2 is due to the fact that the order of the derivative in vr is two less 
than the total order. 

In the scale of weighted spaces 

X" = {{R,7r) I R G An{t){L'^Hl/^ n L^I^^^.e)^^ ^ An{t)L]} (2.243) 

and with the notation 

mi + . . . + rUn = m, (2.244) 
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one also has that 



A„+„_2(i)L2w-l/2,6/5 < C1C2 ||(5-Rl,(57ri)||x'"l • ■•• • \\{SRn,STTn)\\x^n 



(2.245) 



Then, in the Taylor expansion of F{Tr, R), where vr = vr'^ + vr^ + . . . and R = R^ + R^ + . . . 
are functions of Rq, the n-th order term is 



F^■((i^"^7r"l),...,(i^"^7^"0)• 

j=l ni+...+nj=n 



(2.246) 



By writing only the n-th order terms in the Taylor expansion of (|2.207p . we see that i?" 
and vr" fulfill the equation 

n 

idtR"" + n^o{t)R'' = Y Yl i^^'((^"^7^"^),•••,(-R"^vr"0) 

j=2 ni+...+nj=n 



E E (aiKo)(i)(vr"S...,vr"Oi?" 

j=l ni-\-...+nj+n=n 



(2.247) 



/" = E E vr"^), • • • , (i^"^ vr-0). 

j=2 ni+...+nj=n 

with initial conditions i?"(0) = h"-F+{Wo) and 7r"(0) = 0. 

Since the expansion is around zero, note again the cancellations of all terms containing 
i?0 or vrO; 

|(.o,^o)= 0, F} 1(^0,^0)= 0, (a,y,o)(t)7r"i?o = 0. (2.248) 

This means that highest-order terms, depending on i?" and vr"", are completely absent from 
the right-hand side, which is consequently given by combinations of lower-order terms. 
We assume that the induction hypothesis 

(2.249) 



||(i?'-,7r-)||^™-i < A;iA:^||i?o||^, 



/2 



holds for every m < n and some constants ki and k2 to be established later. Then the 
right-hand sides of equations (|2.247|) obey the bounds 



j=2 ni+...+nj=n 



< 



(2.250) 



<Ci^(C72fci)^fc?po||^i 
i=2 



/2 



and likewise for the other two terms, that is 



Y E (5iKo)(t)(vr"S...,vr"0^" 

j=l ni+...+nj+n=n 



< 



(2.251) 



<CiY(C2kiyk^\\Ro\\l,/„ 

j=2 
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respectively 



^ ^ i^■((i^"^7r"l),...,(i^"^7^-0) 

j=2 ni + ...+nj=n 



< 

(2.252) 

<Ci^iC2hyk^\\Ro\\]j,^,. 

J=2 



Importantly, summation starts from j = 2 on the right-hand side, because the terms 
involving first-order differentials vanish, following (|2.248p . 
The powers of t add up in correctly, since 

n 

j -2 + ^{nj -1) = n-2 <n-l. (2.253) 

i=2 

By making ki sufficiently small, we can get the sum X^"=2(^2^i)-' to be bounded and 
small without regard to n and thus we can put 

Cklk^WRoir^,^, (2.254) 

on the right-hand side of (l2.25Ufl2:252]) . 

At this point we solve the equation system (I2.247P for i?" and tt" in the — by now — 
customary manner. For the modulation path tt", (|2.252p is exactly the estimate we need. 

We split i?" into its projections on the three parts of the spectrum, absolutely contin- 
uous, null, and imaginary, and estimate each separately in the weighted Strichartz norm. 

As previously (see (j2. 65112. 68)) ). we first apply a unitary transformation U{t) to the 
equation, such that 

U{t) = e/o + ^((«°)'(«) - \v'is)\^)as) ds 

Z{t) = U{t)R''{t) (2.255) 
Ty(7r°) = U{t)W^o. 

In case of (I2.247p . this transformation takes a particularly simple form, as and are 
constant. 

For the projection on the continuous spectrum Pc{t)Z{t), Strichartz estimates lead 
directly to 

ll^^(^)^(*)IU„W(LrHy^nL?H.i/^-) ^ Ckfk^moW],.;.- (2.256) 
Concerning the projection on the imaginary spectrum, note the convolution estimate 

f\-^('-^)An{s)ds<CAn{t), 

V ^ (2.257) 

poo /"OO V / 

/ e'''^^-''Un{s)ds= / e-''''An{t + s)ds<CAn{t). 
Jt Jo 

Since the right-hand side manifests only polynomial growth, it follows that there exists 
a unique subexponential solution to the ordinary differential equation that describes the 
projection on the imaginary spectrum, corresponding to a suitable value of the parameter 
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Then we obtain 

||P.„.(t)^(t)IU^(,)(,^^i/.,^.^i/.6) < CklmRor^.;.. (2.258) 
The orthogonahty condition is fulfilled approximately, in the sense that 

n 

{R\Ef{W^o))+Y^ (i?",(5^H;(Ty,o))(7r"\...,vr"0) = 0. (2.259) 

j=l 1^1+. ..+nj+h=n 

Under our induction hypothesis, this is enough for an appropriate bound on the Pq com- 
ponent: 

ll^o(*)^(*)llA4t)(Lr=iTy^nL?^i/^'«) ^ Ckjk^WRornr/.. (2.260) 
In conclusion, for the unique suitable value of the parameter /i", 

||(i?",7r")||^„-i < Cklk^WRoWl,,,. (2.261) 

By setting ki sufficiently small, we obtain 

||(i?",7r")||^„-i <A;iA;^||i?o|i:^i/2. (2.262) 

In order to compensate for the smallness of ki, we also need to set /c2 to be large, so that 
the initial condition for induction 

\\iR\7r^)\\xo <kik2\\Ro\\H^/2 (2.263) 

is verified. 

Next, we verify that i?" and vr" are indeed the n-th order terms in the power series 
expansion of R and vr. Denote 

= + i?i + . . . + i?", = 7r° + + . . . + vr". (2.264) 

We obtain that, up to an error O"^^ of order n + 1, S'"' and S" solve an equation system 
of the form (j2.226jl : 

idtS"" + (t)5" = i^(S", 5") + 0"+^ 

(2.265) 

/e.. = F/(S",5") + 0;+\ fe{a,T,Vk,Dk}, 

with error terms of size 

ll^"^'ll(i>-L?VK™^ + \\OT%)^L} < C\\Ror+%. (2.266) 

Following a comparison in the exponentially weighted space Y between (S'"',S") and 
{R,tt), we obtain that 

||(i?,vr) - (5",S")||y < C{n)\\Ro\\f,%. (2.267) 

This concludes the proof of analyticity for R and vr. 

Finally, recall that, by (I2.99p . h{RQ,Wo) has the formula 

hiRo,Wo) = - / e-/o-(^^W)'^-((F,a3F-(W^^(t)))- 

-3d(0(a(t))-^i?,icT3F-(W^^(t)))+ (2.268) 
+ {a{t)r^{R,ia3{d^F~{W^{t)))ir{t))) dt. 

Since all the components that enter this formula are analytic and grow in time more slowly 
than e**"i , we directly obtain a power series expansion for h. 
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Alternatively, we can write h = + + . . . and find each term /i" as the initial data 
of {R^,it'^) in equation (j2.247p . /i" is the unique value that makes the solution bounded. 
By ()2.268p . it follows that for each n /i" is given by an n-linear form in Rq, whose norm 
grows only exponentially with n. 

Following a comparison between (S'",S") and (i?, vr), we obtain as a byproduct an 
estimate for the difference 

I/, _ (/,0 + . . . + /,n)| < c\Cn\\Ro\\f,%. (2.269) 

This gives us another way to establish the analyticity of h. □ 

2.6. The centre-stable manifold. Finally, we state the connection between the mani- 
fold M issued by our proof and the centre-stable manifold of [BaJo]. 

Proposition 2.16. M is a centre-stable manifold in the sense of Bates-Jones. 

Proof. To begin with, we rewrite equation (10. ip to make it fit the framework of the theory 
of Bates- Jones |BAJoj . 

Consider a soliton WTrg{t), described by the constant path 7ro(t) = (1,0,0,0) (without 
loss of generality) . Thus 

w„{t) = (;"fj(:\')) (2.270) 

for all t. Linearizing the equation around this constant path yields, for R = ^ — W-^q, 
equation (j2.17p . which in this case takes the form 

dtR - in^oR = ^(^' W^o)- (2.271) 

Making the substitution 

Z = e-'^^'^R, (2.272) 

we see that 

idtZ + HiWiTro))Z = N{Z, W{Tro)), (2.273) 

where 



and 



Note that the right-hand side terms are at least quadratic in Z, due to linearizing around 
a constant path. 

The spectrum of W is (t('H) = (— oo, — 1]U[1, oo)U{0, ibicr}; for a more detailed discussion 
see Section [2.21 The stable spectrum is — icr, the unstable spectrum is ia, and everything 
else belongs to the centre. 

One can check that all the conditions of jBAJo] are met in regard to ()2.273p , leading to 
the existence of a centre-stable manifold in this setting. This was done by [GJLSj . whose 
main theorem we cited in the introduction. 

Indeed, in the Banach space X = ff*, s > 3/2, Ti is a closed, densely defined operator 
with the required spectral properties. The nonlinearity N{Zi,W{tto)) has the Lipschitz 
property with arbitrarily small constant, since H'^ C L°° is an algebra: 

\\N{Zi,W{7To)) - N{Z2,W{7TomH^ < C max( || Zi ||^^. , || Z2 ) II ^1 " ^2 1| (2.276) 
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We have also exhibited another stable manifold, M C H^/'^ given by Definition 12.21 
which is invariant under symmetry transformations and under the time evolution induced 
by (jO.ip . Our claim, which we prove in the sequel, is that ^' belongs to this manifold if 
and only if Z belongs to the other one, that is 

M = M-^, ^= (2.277) 

is a centre-stable manifold for (j2.273p relative to a neighborhood V of 0, namely V = {Z \ 
\\Z\\iji/2 < 5o} for some small Sq. 

The main difference between the two results is that i/^/^ is not an algebra, so the 
conditions for Bates-Jones's result, Theorem 11.4^ are not met in this space. Indeed, if 
ip S H^/"^, it does not follow that IV'PV' 6 H^^"^ , much less that this nonlinearity is Lipschitz 
continuous (see Condition H3 preceding Theorem 1 1.4p . 

However, even though the hypothesis fails in iJ^/^, the conclusion of [BaJo] — the 
existence of a centre-stable manifold M — still holds, following a proper use of Strichartz 
inequalities. 

To show this, we verify the three properties listed in Definition 11.11 N is t-invariant 
with respect to a neighborhood of 7t^^{N') contains a neighborhood of in X'^ © X^, 
andAAnVF" = {0}. 

The t-invariance of J\f relative to V follows from definition (see the introduction) and 
Proposition 1 2 . 1 3l Indeed, the invariance established by Corollarv l2.14l is strictly stronger 
than i-invariance, as it holds globally in time. 

The fact that 7r'^'^{J\f) contains a neighborhood of in X'^ © X^ is a consequence of the 
local invertibility of J- established in Proposition 12.111 

Finally, we show that MCiW^ = {0}, where is the unstable manifold of the equation 
(the set of solutions that decay exponentially at — oo). To this purpose, in the sequel we 
follow the proof of |Bec1] with slight modifications. 

Consider a solution Z° e iV n VF" of (|2.273p . meaning that Z° is defined for alH < 0, 

||^°(t)||^i/2 <<^o (2.278) 

for some small 5q and all negative t, and Z^ decays exponentially as t —oo, meaning 
that for some constants Ci and C2 and t < 

||^°(t)|lHi/2 <Cie^'* (2.279) 

(even though polynomial decay is sufficient for a contradiction). Using Strichartz estimates 
for the free Schrodinger evolution, we obtain that for every T 

\\^°huT,T+i)w^^^'' - ^ll^°llL-(r,T+i)i/y2- (^-^^O) 

Therefore, for T < 0, 

\\Z'\\LU~oo,TW^^^^^^Ce^''^- (2-281) 

Thus, ^(t) = e^^"'^Z^{t) + Wnoit) is a small asymptotically stable solution of (jO.ip in a 
sense similar to Definition 12.41 but as t goes to —00. 

Therefore, for t < one can write ^ = R{t) + WTrit), such that the orthogonality 
condition 

Po{W^{t))R{t) = (2.282) 
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is satisfied and for T < 

||i?(r)||^,/. < C6, PII,.(_,^]^y^nL?(-oo,T]H'i/^'« ^ Ce^'^- (2-283) 
The path vr satisfies the corresponding modulation equations and 

||^ILi(-oo,T] < Ce^'""- (2-284) 

Since ^(0) G A^, the path vr extends to ah positive times and R extends to a solution 
bounded in the Strichartz norm: 

rTorn N/rl/2^r2rn \Ti^l/2.6 ^ C6. (2.285) 

As in previous arguments, let U{t) be the family of isometries defined for vr 

U{t) = e/o(2^('5)^ + ^("^(•s) - \v{s)?)<yz) ds (2.286) 

and set 

Z{t) = U{t)R{t). (2.287) 

Then Z{t) satisfies the equation 

idtZ -n{t)Z = F. (2.288) 

Decompose Z into its projections on the continuous, imaginary, and zero spectrum of 
n{t) and let 

^(.T) = \\Z\\l^(^_^^t]h1/2^lU-oo,t]w^^'''^ + MlI{-<^,t]- (2-289) 
Observe that 5{t) ^ as t ^ — oo, so we can assume it to be arbitrarily small. 
By means of Strichartz estimates one obtains that 



<C5(T)||^||,.(_^^j^y2, 



(2.290) 



because the right-hand side contains only quadratic or higher degree terms. 
Let 

P,^Z{t) = b_{t)F"{t)+b+{t)F+{t). (2.291) 
Since Z{t) is bounded as t — > — oo, we can use Lemma 12.41 for t — > — oo and obtain 

(2.292) 



h+{t) = e(*'^)"6+(T) - r e^'-'^''N+{s)ds. 



Therefore 



WP^-^hn-oo^TWl'^ ^ C{\\P+Z{T)\\^.j, + 5(T)||Z||^^^(_^^j^v2)- (2.293) 
We have constructed Z such that the orthogonality condition holds, so 

Po{t)Z{t) = 0. (2.294) 
Putting these estimates together, we arrive at 

\\^\\L^i-oo,T]Hl^^ ^ ^(^(^)ll^llL?o(-oo,T]i^y^ + \\P^Z {T)\\ ^.,,) . (2.295) 

For sufficiently negative Tq, it follows that ||Z(t)||^i/2 < C\\P+Z{t)\\^i/2, for any t < Tq. 
The converse is obviously true, so the two norms are comparable. 
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Furthermore, by reiteration one has that 

11(1 - < C6it)\\P+Z{t)\\^,;,. (2.296) 

Next, assume that Z{0) is on M, meaning that Z(0) + $ G M. 

If the size 5o that appears in Definition 12.41 of M is sufficiently smah, it follows that 
II Z(t) 11^1/2 is bounded from below as t — > cxo. Indeed, to a first order, Z{t) is given by the 
time evolution of the initial data, following (j2.218p . and this expression cannot vanish. 

On the other hand. Lemma 12.41 implies that 

eit-^)-\N+{s)\ds (2.297) 

and thus ||P+Z(i)||^i/2 goes to zero and can be made arbitrarily small as t ^ oo. 

Lemma 2.4 of [BaJo| implies, under even more general conditions, that if the ratio 
||P+Z(To) 11^:^1/2/11 (/ — P+)Z(Tq)\\^i/2 is small enough, it will stay bounded for all t < Tq. 
The proof of this result is based on Gronwall's inequality. 

However, this contradicts our previous conclusion (j2.296p . stating that 

11(1 - P+)Z(t)||^i/2/||P+Z(t)||^v2 < C6it) (2.298) 

goes to as i goes to —00. Therefore, Z can only be 0. 

This proves that M D = {0}. In other words, there are no exponentially unstable 
solutions in J\f in the sense of |BaJo] . The final requirement for to be a centre-stable 
manifold in the sense of Definition II . II is thus met. □ 

2.7. Scattering. Strichartz space bounds and estimates imply that the radiation term 
scatters like the solution of the free equation, meaning 

r{t) = e-''^rfree + 0^1/2(1) (2.299) 

for some r j^ee 

As a reminder, R satisfies the equation (12.17P 

idtR-n^{t)R = F, (2.300) 

where TiT^^t) = Ao"3 + Vn-(t), and R has finite Strichartz norm, 

while F has finite dual Strichartz norm, 

||P||^2vj/i/2,6/5 < 00. (2.302) 

Rewrite (|2.300l) as 

idtR- Aa3R = F -V^it)R. (2.303) 

By Duhamel's formula, 

R(t) = e-'*^"^P(0) - i f e-^(*-^)^"^ {F{s) - V^{s)R{s)) ds. (2.304) 

JO 

/■oo 

Rfree = R{0)-i e'^^"'' {F{t) - V^{t)R{t)) dt. (2.305) 
Jo 

/oo 
gisA<x3 _ dt. (2.306) 



Let 



Then 
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Note that 

||F - r2ii/l/2,6/5 < 

11 11 11 11 11 11 (2-307) 

imphes 

^hm \\x[t,oo){s)iF{s) - V^is)R{s))\\^,^if,,e/B = 0. (2.308) 

g-iiA<73 i^ging isometry, it follows that R{t) - e"'^^"^ Rfree ^ in H^/'^. 

This leads to the desired conclusion, upon passing to the scalar functions r and rj^eei 



where R = [_ \ and i?/ree 

. ^ / V f free 



ree 



3. Linear estimates 

We seek a dispersive estimate for the linear time-dependent equation. The highest-order 
terms we need to take into account are of the form v{t)VZ{t), where v{t) is small in the 
norm. Z is a solution of the Schrodinger equation and thus has a finite Strichartz norm. 

3.1. Notations. Consider the linear Schrodinger equation in M.^ 

idtZ + nZ = F, Z(0) given, (3.1) 

where 

Wi and W2 are taken to be real-valued. We also assume that Wi and W2 possess a 
sufficiently large amount of polynomial decay and smoothness, meaning that {x)~'^d'^V G 
for all m, n < N and some sufficiently large N. 

All subsequent computations take place under these far from optimal conditions. A 
derivation of similar results in a sharp (scaling-invariant) or almost sharp setting can be 
found in [Bec2| and |Bec3 ] . 

In keeping with the nonlinear problem, we assume here as well that TL has no eigenvalues 
or resonances in (—00, —fj] U [fi, 00). 

For simplicity, the entire subsequent discussion revolves around the case of three spatial 
dimensions. We focus on H^^^ estimates, involving half-derivatives, which are of use in 
the study of the the nonlinear problem under consideration. 

Following [Sch] and [ERScj . we list the spectral properties of such operators. 

Proposition 3.1. Let H = Hq + V be as in 13. 2|) and assume that {x)~^d'^V G L°° 
for all m, n < N and some sufficiently large N . Then o'{TC) C M U and o'ad'H) = 
(-00, -/u] U [/i,oo). 

Assume that TC has no eigenvalues or resonances in (— co,— /u] U [/U,oo). Then point 
spectrum is simple, with the possible exception of the eigenvalue at zero. The finitely 
many Riesz projections P^. , ^ < j < n, corresponding to the point spectrum are given by 
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have finite rank, and their ranges, as well as those of their adjoints P^, , are spanned by 
Schwartz functions. The continuous spectrum projection Pc is given by 

Pc = I-Y.PCr (3-4) 

i=i 

We also state the limiting absorption principle for 7i as follows. 

Proposition 3.2. Let Ti = T-Lq + V he as in \'S.2\) and assume that {x)^"^d^V G L°° for 
all m, n < N and some sufficiently large N . Assume that there are no eigenvalues or 
resonances ofTi embedded in {—oo,—fi] U [fi,oo). Then 

sup ||i?y(A ± i0)||^6/5_i6 < oo, sup ||i?y(A ± i0)\\^i/2,6/5^^ri/2,6 < oo. (3.5) 



Proof. Write V = V1V2, where 
Since 

Ry{\) = i?o(A) - i?o(A)Fi?o(A) + Rg{\)Vi{I + V2Ro{\)Vi)-^V2Ro{\). (3.7) 

the boundedness of iiy(A it iO) follows from the existence of a bounded inverse for I + 
V2-Ro(A±iO)yi as an operator from L? to itself (or iiT^/^ to itself). By Fredholm's Theorem, 
this is implied by the nonexistence of a function f £ (or H^/'^) such that 

/ + y2i2o(A±iO)T/i/ = 0. (3.8) 

Indeed, Vif would have to be an eigenfunction or resonance for at A € (—00, — /i]U[/i, 00), 
contradicting our hypothesis. Note that only it/i could actually be resonances. 
Finally, the uniform boundedness of the norm follows from the fact that 

lim ||y2^o(A±iO)yi||i2^i2 = (3.9) 

A^oo 

and likewise in H^l'^. □ 
The resolvent of the unperturbed Hamiltonian, i?o(A) = {Hq — A)^^, has the kernel 

i?o(A^ + /^)(x,y) = - . (3.10) 

i?o(A) = (Hq— A) ^-"^ is an analytic function on C\((— 00, —/^]U[/^, cxd)). It can be extended 
to a continuous function in the closed lower half-plane or the closed upper half-plane, but 
not both at once, due to a jump discontinuity on the real line. 

3.2. Strichartz estimates. We begin by recalling the following basic connection between 
the free evolution e**^" and the resolvent. 

Lemma 3.3. 

l|e**^1lLi£{LinL2,L-+L2) < 00. (3.11) 
Furthermore, the Fourier transform of e**^" is well-defined, in the sense that the integral 

lim r e-**^e**^o dt (3.12) 
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converges, for A in the lower half-plane, in the C{L^ n L^,L^ + operator norm to 
iRo{X) and to iRo{X — iO) for A G M. For real A 

lim [\-'^^e'*'^'>dt = i{Ro{X-iO)- Ro{X + iO)) (3.13) 

also converges in the operator norm. 

Proof. To show (13.1ip . simply use the dispersive estimate 

||e^*^o||ii^ioo < Cr3/2 (3.14) 

for t > 1 and the unitarity of the evolution 

||e**^o||^2^i2 <C (3.15) 

for t < 1. Then the integral in (I3.12p is well-defined and converges in the operator norm 
due to dominated convergence. Furthermore, 

r e-''^^-''h''^^ f dt = iR^iX - ie){I - e-^p(^-*^)e^^^«)/. (3.16) 
Jo 

Since / G L^, the right-hand side converges to iRo{X — ie)f . Letting e go to zero, the 
result also follows for A on the real line, by dominated convergence. One must use the 
explicit form of the kernel (jS.lOp . 

Finally, (j3.13p follows by separately handling the integrals from — p to and from 
to p. □ 

This admits the following generalization to the case where a potential is present. 

Lemma 3.4. Assume V G and the Hamiltonian 7i = Tio + V is described by ( 13. 2|) . 
Then the equation 

idtZ + nZ = F, Z(0) given, (3.17) 
admits a weak solution Z for Z(0) G L^, F G L^L'^ and for t > 

||^(t)||2<Ce*ll^ll-||Z(0)||2+ re(*-^)ll^ll°°||F(s)||2<is. (3.18) 



Furthermore, Rv{X) is the Fourier transform o/e**^.' forlmX < —Halloo 

lim r e-'^^e'^'^ dt = iRv{X). (3.19) 

Proof. We introduce an auxiliary variable and write 

idtZ + T-CqZ = F - VZi, Z(0) given. (3.20) 

Over a sufficiently small time interval [T, T + e], whose size e only depends on ||l^||oo5 the 
map that associates Z to some given Zi is a contraction, in a sufficiently large ball in 
L'^L'^. The fixed point of this contraction mapping is then a solution to ()3.17p . 

This shows that the equation is locally solvable and, by bootstrapping, since the length 
of the interval is independent of the size of F and of the initial data Z(T), we obtain an 
exponentially growing global solution. The bound (j3.18p follows by Gronwall's inequality. 

(fSTTOl) follows in the same manner as (j3. 12p . □ 

In order to obtain Strichartz estimates in the time-independent case, we first derive the 
following representation formula (also proved in |Sch| under more restrictive assumptions). 
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Lemma 3.5. Let V £ H L°° and Ti be given by {3 .21) . Assume that there are no 
eigenvalues of H embedded in its essential spectrum (— oo, — //] U [/U, cxd) and only finitely 
many outside of it. Then for f , g G L"^ 

{Pcf, 9) = !^ I {{Rv{X- iO) -Rv{X + iO))f, g) dX, (3.21) 

ZTTl J(_oo,-/i]U[At,oo) 

in the sense that the expression under the integral is absolutely integrable and its integral 
equals the left-hand side. 

Proof. Let V = V1V2, where 

^^ = ^H^2 '^2=Ih^2 '^3=(o _J. (3.22) 



For y > \\V\\oc, I + y2Ro{X±iy)Vi must be invertible. Indeed, Vi and V2 are bounded 

1 f e~v^l^l 



operators of norm at most and 



1 r e~vym 

\\Ro{X±iy)\\2^2<^ -^^dx = l/y. (3.23) 

4vr Jm.3 \x\ 

Therefore one can exphcitly construct the inverse (/ + V2Ro{X it iy)Vi)~^ by means of a 
power series. 

By Lemma 13.41 Y^>n(e^*^e~^* f. a) is an exponentially decaying function of t and its 
Fourier transform is 



(e-(^+*")*e**^/,5) dy = -i{Rv{X - iy)f,g). (3.24) 

Combining this with the analogous result for the positive side, we see that 

{{e'''^e-y\'\f,g)r = i{{Rv{X + iy) - Rv{X - iy))f,g). (3.25) 
The right-hand side is absolutely integrable, because 

Rv{X) = Ro{X) - Ro{X)VRo{X) + Ro{X)Vi{I + V2Ro{X)Vi)-W2Ro{X) (3.26) 

and 

OD 

\{{Ro{X + iy)-Ro{X-iy))f,g)\dX 

■oo 

/oo 1 
-(((i?o(A + iy) - Ro{X - iy))f, f) + {{Rq{X + iy) - Ro{X - iy))g,g)) dX (3.27) 



2 

The remaining terms are absolutely integrable due to smoothing estimates: 

\\\V\"^RQ{X±iy)f\\ldX<C\\f\\l (3.28) 
By the Fourier inversion formula, (j3.25p implies 

/ x{X){{Rv{X + iy)- Rv{X-iy))f,g)dX = {f,g). (3.29) 

^TT ./in. 
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We then shift the integration contour toward the essential spectrum (— oo,— /u] U [//,oo), 
leaving behind circular contours around the finitely many (by Fredholm's Theorem) eigen- 
values. Each contour integral becomes a corresponding Riesz projection. 

What is left is the projection on the continuous spectrum. The integral is still 
absolutely convergent due to (|3.26p . (j3.27p . and smoothing estimates. (|3.2ip follows. □ 

Using this representation, we proceed to derive H^/"^ Strichartz estimates in the time- 
independent case. 

Theorem 3.6. Let Z he a solution of the linear Schrddinger equation 

idtZ + nZ = F, Z(0) given. (3.30) 

Assume that 7i = TCq + V, {x)^d^V G for all m, n < N and sufficiently large N, 
that V is matrix-valued as in i t3.2|) . and that no eigenvalues or resonances are present in 
(— oo, — /u] U [/i, cxd). Then 



W^'^^Wh^Hl'^nL^^wl'^''' - <^(II^(0)IIhi/2 + ll^llLii/i/2+L2iy™5j- (3-31) 

Proof. Let F € L\{L]j. n L^), G G L'^{L\ n L^) have compact support in t and consider 
the forward time evolution 

{TF){t)= f e'^^-'^'^PcF{s)ds. (3.32) 

Jt>s 

TF(t) is in for all t and grows at most exponentially, so its Fourier transform is well- 
defined for ImA < — ||F||oo (where, in particular, ||i?y(A)||2^2 is bounded): 

TF{X) = iRv{X)PcF{X). (3.33) 
Using the representation formula p.2ip for Pc given by Lemma [331 for f^ 9 ^ 

{Rv{Xo)Pcf,g) = ^ [ {Rv{Xo){Rv{X-tO)-Rv{X + iO))f,g)dX 

= -^/ /^—(Ry(X-iO)-Rv{X + iO))f,g)dX. 

^^W(-CXD,-At]U[/i,00) ^-^ - -^0 ' 

Here we used the resolvent identity: for all Ai, A2 in the resolvent set, 

RviXi) - Rv(X2) = (Ai - X2)Rv{Xi)Rv{X2). 
From the absolute convergence of the integral for f,g£L^,it immediately follows that 

\\Rv{Xo)Pc\\l^^l^ < -T7T—, TV- (3-35) 

d(Ao, (-CX), -^J U [^,00)) 

The limiting absorption principle provides that 

sup \\Ry{X ±iQ)\\]^e,/h ^1^6 < 00 (3.36) 

Ae (— 00,— ^]U[/^,oo) 

and thus, for 6/5 < p < 2, ({Rv{X - iO) - Rv{,X + iO))f,g'^ G ^2p/(5p-6)_ Therefore, for 
Ao {—oo, —fj] U [fi, 00) and 6/5 < p < 2, 

\\Rv{Xo)Pc\\lp^lp' < (3-37) 

Also note that 

RyPc = RqPc - RqVPcRo + RqVPcRvVRo (3.38) 
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and consequently, for Aq (— 00, — /u] U cxd), 

\\Rv{^0)Pc\\L(i/5-,L(i < ||-Ro(Ao)|Il6/5^2.6 + l|-^o(Ao)||^6/5^i6ll^llL3/2 + 

II^o(Ao)||li^l6 11^1123/2 Py(Ao)^'c|| 

(3.39) 

Using the limiting absorption principle, it follows that 

sup \\Rv{X)Pc\\l^/5^l(^ < 00. (3.40) 
Aec 

For y > \\V\\oo, both e-y\TF){t) and e?'*G(t) belong to Ll^. Taking the Fourier 
transform in by Plancherel's theorem 

{(TF)it),G{t))dt = ^ I {{e-y\TF){t))\{ey'G{t))'')d\ 

■'^ (3.41) 
= 7^ / {Rv{\-iy)PcF{\-iy),G^){\-iy))d\. 

Here (•, •) is the real dot product. 

The pairing makes sense because RyPc is analytic in the lower half-plane, as a family 
of bounded operators from L^/^ to . Furthermore, F and G are analytic due to the 
compactness of the supports and, for every y G M, F{\ + iy) € L\L^J^ and likewise for G. 

We shift the integration line toward the real axis, obtaining 

/ {{TF){t),G{t))dt = ^ [ {RviX-iO)PcFiX),G(Pt){\))d\. (3.42) 



This leads to 



/ {{TF){t), G{t)) dt < Gsnp \\Rv{X - iO)Pc\\aL^/, l^>) 11^11^2^6/5 \\G\\ , e/. 



(3.43) 



Then, we remove our previous assumption that F and G should have compact support. 
This establishes the Strichartz estimate 



eiii-^)'Hp^F{s) ds 

t>s 



<C||F|| 6/5. (3.44) 



Using Duhamel's formula we obtain all the other Strichartz estimates. 

For H^^^ Strichartz estimates, we start again from (j3.42p and use the fact that 

sup \\Rv{X - iO)Pc||£(vi/i/2,6/5^iyi/2.6) < oo. (3.45) 



Its derivation is exactly analogous to that of (j3.40p . □ 

3.3. The time-dependent case. We turn to time-dependent equations. We prove a 
concrete result in the case of interest, while at the same time placing it within a more 
general framework. 

Given parameters A(t) and v{t) = {vi{t),V2{t),V3(t)), consider the family of isometrics 

U{t) = e/o(^'^(^)^ + «^(s)f^3) ds^ ^3 
The rate of change of U{t) is then controlled by the norm 

\\AmLr + \\vmLr- (3.47) 
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The linearized Schrodinger equation under study has the form 

idtR{t) + {Ho + U{t)-^VU{t))R{t) = F{t), R{0) given. (3.48) 

Observe that the Hamiltonian at time t is U{t)~^(7io + V)U{t), since Ti and U commute, 
i.e. it is Tio + V = 7i conjugated by U{t). 

Let Z{t) = U{t)R{t). We rewrite the equation in the new variable Z, obtaining 

idtZ{t) - idtU{t)U{t)-^Z{t) + noZ{t) + VZ{t) = U{t)F{t), Z(0) = ii(0). (3.49) 

In the subsequent lemma we encapsulate all the properties of U that we actually use in 
our study of (ITiHD and (fXIHD . 

Lemma 3.7. Let U{t) be defined by < I3.46|) . 

(1) U{t) is a strongly continuous family ofW^'^ isometries, for s G M, 1 < p < oo. 

(2) For every t and r > 0, U{t) and U (r) commute with TCq and each other. 

(3) For some N and all < s < 1, there exists e{N) > such that 

|K.)-^([/(t)C/(r)-ie^(*-)«o _ e^(*-^)^«)(x)-^||^(^.^.,^.^.^ < 

<C{\\A{t)U^ + \\v{t)U^)<^\ 

Proof. It is easy to verify the first two properties directly. In particular, the Laplace 
operator A commutes with translations. 
Concerning the last property, we compare 

r(t,s) = e^(*-^)^o (3.51) 

and 

f{t,s) = e«(*-^)WogX,'(2i'(T)V+jA(T)(T3)rfr_ (-3^52) 

On one hand, from dispersive estimates, 

\\T{t, s) - t{t, s)||i^oo < C\t - s|-3/2. (3.53) 

On the other hand, 

\\T{t,s)-f{t,s)\\2-^2<C. (3.54) 

It follows that for iV > 1 

IKx)-^(r-f)(x)-^||^(^.^,^.^) < a (3.55) 

This holds with a constant independent of A and v. 

Finally, we prove the last stated property. Assume A and v are small. We first consider 
the case when there is no translation movement due to v and we only have to handle 
oscillation, due to A. Denote this kernel by Tqsc, that is 

foscit, s) = e*(*~^)^Ogr (iA(r)-3) dr_ (3,56) 

One has 

e*'' - 1 < Cmin(l,a) (3.57) 

and thus 

\\Tit,s) - foscit, s)\\2^2 < Cmin(l, \\A\Ut - s\) < C\\A\Ut - s|^ (3.58) 
Therefore, for large enough N, it follows from ()3.53p and from ()3.58p that 



{X 



■-^(T - r,,e)(x)-^||^(^2 .2 . < C\\A\\t. (3.59) 
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Next, we consider the case when v is not necessarily zero. Let d{t) = £v{T)dT. Then 

(3.60) 



^-i{t-s)A^{J*2v{T)dT)V ^ 



o/a ..Y l^-vr (x-y)(d(t)-d(s)) (d(t)-d(s)y \ 



Ax-y)(d(t)-d(s)) 

e ~- - 1 



(-47ri)3/ 

. (d(t)-d(s))^ 

We treat the last factor *-« in the same manner in which we treated the factors 

. (d(t)-d(»))^ 

containing A. Consider the kernel that contains those factors together with e , 

(x-y)(d(t)-d(s)) 

leaving aside e : 

fi{t, s) = e^ft-^)Hoe/; iiAiT)a3) dT /'"-\_f^^^ _ (3.61) 
For this kernel as well, one has 

||(x)-^(r-ri)(x)-^||^(^.^,^.^) < C{\\A\\t + ll^^ll^^). (3.62) 
Considering the fact that 

' < Cmin(l, ||i;||oo(N + \y\)) < C\\v\\U\x\ + \y\r, (3.63) 
it follows that for large enough N 

\\{x)-''ifit,s)-m,s)){x)-^\\2^2 < C\\v\Ut-s\-'/''. (3.64) 
We also have the trivial bound 

\\f{t,s)-fi{t,s)\\2^2<C. (3.65) 

Therefore 

||(x)-^(f -fi)(x)-^||^(^2^,^2^) < C\\v\\t. (3.66) 

Overall, we find that 

|Kx)-^(r - r)(x)-^||^(^.^,^.^) < C7(||^||^=^ + \\v\\t). (3.67) 
Finally, note that 



|9fc(x)-^(r - r)(x)-^F|L.^ < c(|Kx)-^-i(f - T){xy'^F\\^. + 



(3.68) 



+ IKx)-^^(T-r)(x)-^^-^F||^.^+ 

+ ||(x)-^(T-r)(a;)-^5,F|L.J 
<C7(||^||^3 + ||^||^^)||F||^.^,, 

in view of the fact that C L^. Consequently, for large enough N, 

||(a;)-^(T - T){x)-^F\\^,^, < C{\\A\\i' + ||^||^3)||F||^.^,. (3.69) 

Interpolating between this and L^, we obtain boundedness on L^H^, < s < 1. □ 



Before proceeding with the proof of the Strichartz estimates in the time-dependent case, 
we require the following technical lemma. 
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Lemma 3.8. Consider H = Hq + V as in < [3 . 2|) such that {x)^d"^V G for every m, 
n < N and some sufficiently large N . Assume that 7i has no eigenvalues or resonances 
embedded in a{7io). Then there exists a decomposition 

V-Ppn = F^{V)F2{V), (3.70) 

where Pp = I-Pc, such that Fi{V) G C{L^, {x)-^L'^), F2{V) G C{{x)^ , L'^) . Moreover, 
it is the case that 

Fi{V) G C{H^/\ {x)-^H^/^), F2{V) G C{{xfH^/^, H^''^). (3.71) 
Proof. Simply take for some large n 

F,{v) = {V - Ppn){xr 

This works because both the range of PpTi and that of its adjoint are spanned by Schwartz 
functions. □ 

Theorem 3.9. Consider equation i f3.48|) . for H = Ho + V as in 113. 2\) : 

idtZ - iv{t)VZ + A{t)a^Z + nZ = F, Z(0) given, (3.73) 

Assume that {x)'^d'^V G L°°, for all m, n < N and sufficiently large N, that ||vl||oo 
and \\v\\oo are sufficiently small (in a manner that depends on V), and that there are no 
eigenvalues or resonances ofTi in {—oo,—fj] U [fi,oo). Then 



WPcZh^LinL^Lt^^ < c[\\mh + m,.Li^,.,r-^) (3.75) 

and 

^^^'^^hrH'J^n^W^/^'" - <^(|I^(0)IIh1/2 + ll^llLli/y2+L2^^1/2,6/5). (3.76) 

Proof. To begin with, the fact that V is regular and of rapid decay guarantees the existence 
of a solution R, albeit one that may grow exponentially. 
As in the time-independent case, let 

Z = P,Z, F = P,F - 2iv{t)[Pc,V]Z + A{t)[Pc,a3]Z. (3.77) 

The equation becomes 

idtZ - iv{t)VZ + A{t)a3Z + nZ = F, Z(0) = Pc^(O) given. (3.78) 

The commutation terms 

2iv{t)[P„V]Z, A{t)[P„a3]Z (3.79) 

are small in the dual Strichartz norm (for small ||f ||oo and ||^||oo) and thus can be con- 
trolled by Strichartz inequalities and a simple fixed point argument. 
Lemma 13.81 then provides a decomposition 

V-Ppn = Fi{V)F2{V), (3.80) 

where Fi(y) and F2{V) are bounded from to {x)~^ L'^. 
Denote 

^F.(y),Fi(y)i^(t) = f F2{V)P,e'^'-'^'^'U{t)U{s)-'F,{V)F{s)ds, (3.81) 
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respectively 

fp^iy),iF{t)= f F2{V)P,e'^'-'^'^^^U{t)U{s)-'F{s)ds. (3.82) 

By Duhamel's formula, 

F2{V)Z{t) = ifF^^v),F,{v)F2{V)Z{t) + fF^^y)j{-iF{s) + <5,=o^(0)). (3.83) 
We compare the time-dependent kernel Tp^(y^^p^(y^ with the time-independent one 

TF,(v),F,iV)F{t) = f F2{V)P,e'^'-'^'^^F,{V)F{s)ds. (3.84) 
By Lemma 13.71 we obtain that 

||A||^-.o - 2>2(v),Fi(y)ll£(L2,,L2 J = 0. (3.85) 

The operator / — ^7Fi(V'),F2(y) is invertible in C{L^^^, L'f .j.). Indeed, its inverse is nothing 
but 

(I - iTF^(^v),F,iV))~'F{t) = F{t) - i f F2(y)P,e^(*-^)^Fi(y)F(s) ds. (3.86) 

The right-hand side is in C{L1^,L1^) due to the Strichartz estimates of Theorem 13.61 
The Duhamel formula proves that it is the inverse of the left-hand side, as claimed. 

It follows that, when ||^||oo and ||u||oo are small enough, I —iTp^^y■^^p^(y^ is also invertible 
in C{Lf ^, This immediately implies the desired Strichartz estimates, as 

Z = {fjj + fp^iy)j{I - iTp^^v),F,iV)r'f,^p^^y)){6t=oZ{0) - iF). (3.87) 

In regard to the iiT^/^ case, we have the decomposition 

V -Pp'H = Fi{V)F2{V), (3.88) 

where 

Fi{V) e £(iJi/2^ {xy^H^I'^), F2{V) e C{{xfH^I'^, H^/^) (3.89) 
respectively. By Lemma 13.81 it follows in the same manner that 

l|Al|r-0 - Tf2{V),F^{V) WciLlHl'^LlHl/^) = (^"^^^ 

Therefore we can invert rp2(y),Fi(y) in the operator algebra C{L^Hx^'^ , L^Hx^'^) and the 
proof of i/^/^ Strichartz estimates proceeds by (|3.87p . □ 

Appendix A. Spaces of functions 

Our computations take place mostly in Lebesgue and Sobolev spaces of functions defined 
on M^^^. This corresponds to three spatial dimensions and one extra dimension that 
accounts for time. 

To begin with, we only consider measurable complex-valued functions or ones that take 
values within a finite-dimensional Banach space. In general, dealing with more general 
Banach space- valued functions, one must distinguish between weak and strong measura- 
bility. 
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We denote the Lebesgue norm of / by by ||/||p. The Sobolev spaces of integral order 
p^n,P g^pg then defined by 

^\a\<n ^ 

for 1 < p < oo and 

ll/llw^n.oo = sup ||a"/||oo. (A.2) 

|a|<n 

whenp = oo. In addition, we consider Sobolev spaces of fractional order, both homogenous 
and inhomogenous: 



k-- = ll(V)7llp, respectively ||/||^,, = |||V|/||p. (A.3) 

Here (V)'' and |V|* denote Fourier multipliers — multiplication on the Fourier side by 

(^)« = (1 + |^|2)«/2 and ICI' respectively. 

When p = 2, the alternate notation = W^'"^ or = W^''^ is customary. 
Strichartz estimates may involve mixed space-time norms of the form 

i/p 



f \\f{x,t)r^.,dtj . (A.4) 
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